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ABSTRACT 


Linearized  three-dimensional  lifting  surface  theory 
is  applied  to  a  variety  of  problems  of  Interest  in  hydrofoil 
application.  The  resulting  Integral  equations  are  solved  by 
numerical  techniques  on  a  high-speed  digital  computer  to 
predict  steady  and  oscillatory  loads. 

The  following  cases  are  discussed:  thin  non-planar 
hydrofoils  at  infinite  Froude  number,  the  effect  of  gravity 
waves  on  the  forces  on  an  oscillating  hydrofoil  near  the 
free  surface,  supercavltatlng  hydrofoils  of  finite  span  in 
steady  and  osclllatoi^  motion.  Nmnerlcal  results  for  lift 
and  moment  are  presented  and  comparison  with  experiment 
is  made  whenever  possible. 

Results  of  these  calculations  indicate  that  this 
numerical  solution  is  a  very  practical  and  efficient  way 
to  obtain  accurate  prediction  of  unsteady  loads  for 
flutter  calculations. 
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CHAPTER  I 

LIFTING  SURFACE  THEORy  FOR 
NONFLANAR  SURFACES 

1.1  Introduction 

Before  the  avalllblllty  of  high  speed  digital  computors, 
methods  based  on  linearized  theory  for  calculating  loads  on 
three  dimensional  lifting  surfaces  in  incompressible  fluids 
were  restricted  to  special  plan  form  shapes  which  could  be 
handled  analytically  or  to  approximate  methods  suitable  for 
hand  calculation. 

In  steady  flow,  a  lifting  surface  is  represented  by  a 
vortex  sheet,  or  by  the  horseshoe  vortices  which  make  up  the 
vortex  sheet.  The  boundary  value  problem  is  formulated 
using  the  velocity  potential.  By  the  principle  of  super¬ 
position,  the  upwash  velocity  on  the  wing  is  written  as  an 
integral  of  the  upwash  velocities  from  each  element  of 
the  vortex  sheet.  This  boundary  value  problem,  together 
with  auxiliary  conditions  such  as  the  Kutta  condition,  is 
solved  for  the  distribution  of  vortex  strength  and  thus 
the  distribution  of  lift  on  the  surface. 

Unfortunately  this  Integral  equation  cannot  be  solved 
analytically  for  wings  of  arbitrary  planform.  For  two 
limiting  cases,  wings  of  large  aspect  ratio  and  wings  of 
small  aspect  ratio,  analytical  results  have  been  obtained. 
For  wings  of  arbitrary  planform,  approximate  techniques 
have  been  used.  For  example,  the  wing  is  represented  by 
a  finite  number  of  horseshoe  vortices  located  on  the  surface 
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of  the  wing;  their  strength  is  determined  by  satisfying 
the  known  upwash  condition  at  control  points  on  the  wing. 

A  discussion  of  some  of  these  techniques,  as  well  as 
numerical  and  experimental  results,  is  given  by  Thwaltes, 
Ref.  (1). 

For  planar*  lifting  surfaces  in  oscillatory  motion 
in  an  infinite  uncompressible  fluid,  analytical  solutions 
have  been  given  for  wings  of  infinite  aspect  ratio 
(Ref.  (2)),  and  for  wings  of  circular  planform,  (Ref.  (5), 
(4)).  Some  of  the  approximate  techniques  for  finite  wings 
described  in  Ref.  (1),  such  as  those  of  Palkner  and 
Multhopp,  have  been  modified  for  oscillating  wings.  A 
discussion  of  the  numerical  techniques  for  oscillating 
finite  wings  developed  before  the  age  of  high  speed  comput¬ 
ing  machines  is  given  in  Ref.  (?)•  These  methods  were 
always  a  compromise  between  "accuracy  of  results  and 
computing  labor  required".  (Ref.  (7)) 

Taking  advantaige  of  high-speed  stored  program 
computing  equipment,  Watkins  and  his  associates  ^Ref.  (6), 
(7^  have  developed  a  most  satisfactory  numerical  method  for 
calculation  of  the  lift  distribution  on  a  finite  planar  wing 
in  steady  or  oscillatory  motion  in  an  infinite  subsonic 
flow.  This  me,iiod  makes  a  direct  attack  on  the  Integral 
equation  relating  upwash  to  lift  distribution  on  the 
wing  by  assuming  a  set  of  modes  for  the  lift  distribution 
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with  unknown  coefficients.  The  coefficients  of  these  lift 
modes  are  found  by  satisfying  the  boundary  condition  of  a 
known  upwash  at  suitably  selected  control  points  on  the 
surface.  The  Integrations  required  for  this  method  would 
be  completely  Impractical  to  do  analytically  but  are  quite 
straightforward  numerically. 

Much  of  the  present  work  Is  an  extension  of  these 
Ideas  to  problems  of  Interest  In  hydrofoil  applications. 
Governing  Integral  equations  are  derived  for  the  following 
cases;  In  Chapter  1,  three-dimensional  nonplanar  oscillating 
lifting  surfaces  in  an  Infinite  fluid  or  beneath  a  free 
surface  at  Infinite  Proude  number;  In  Chapter  2,  two-  and 
three-dimensional  oscillating  foils  beneath  a  free  surface 
Including  the  effects  of  gravity  waves  generated  by  the 
motion,  (finite  Proude  number);  In  Chapter  3#  three- 
dimensional  steady  amd  oscillating  supercavltatlng  hydro¬ 
foils  in  an  infinite  fluid.  Numerical  results  are  presented 
for  most  of  these  cases  and  comparisons  are  made  with 
experiments  where  possible. 

These  Integral  equations  are  derived  by  taking  advantage 
of  the  fact  that  In  linearized  theory  for  an  Incompressible 
fluid,  the  perturbation  pressure  satisfies  Laplace's 
equation.  As  an  Introduction  to  these  techniques  we 


shall  discuss  the  Integral  equation  for  a  nonplanar  lifting 
surface  In  an  Infinite  fluid. 


1.2  The  Integral  Equation  for  a  Nonplanar  Lifting  Surface 


In  the  linearized  theory  for  the  flow  of  an  Incompress¬ 
ible  fluid  with  free  stream  velocity  If  In  the  positive 
direction  about  a  thin  lifting  surface  In  steady  or  unsteady 
motion,  the  perturbation  pressure  Is  given  by  a  linear 
operation  on  the  perturbation  velocity  potential,  namely  by 


A 

-E 

p 


-  j-  ^  ^  a/  /, 

We  consider  simple  harmonic  motion  with  frequency  cd 


P  e 


& 


(1.2) 


(1-3) 


If  we  nondlmenslonallze  all  physical  quantities  with  respect 
to  free  stream  velocity  LT  ,  density  ^  ,  and  root 

semichord  h  g  ,  the  complex  amplitudes  of  p  and  ^ 
are  related  by  equation  (1.4) 

-  P  -  1^ 

where  p,  ^  are  now  dimensionless  variables,  k 

Is  the  reduced  frequency 

Equation  (1.3)  Is  solved  for  as  a  function  of 

p  with  the  boundary  condition  oo)  ^  O 
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t'fxj  = 


e 


c  y) 
d  X 


-  (TO 


(1-5) 


In  the  literature  it  Is  common  to  define  an  acceleration 
potential  ^  . 


-  i  Lf  4.  (1.6) 

‘■4^  ^  -  P  (1.7) 


For  convenience,  however,  we  will  work  with  pressure 
directly. 

In  the  linearized  theory  for  an  incompressible 
Irrotatlonal  flow,  both  p  and  If  are  potential  functions 
satisfying  Laplace's  equation. 


VV  =0 

(1.8) 

vV  =  0 

(1.9) 

Qreen's  theorem  (for  a  potential  function)  Is  written 
for  a  closed  surface  surrounding  the  fluid,  Ir  this  case, 
a  closed  Inner  boundary  plus  the  surface  at  Infinity  which 
does  not  contribute  to  the  Integrals  (See  reference  (5)  for 
a  discussion  of  Qreen's  Theorem).  In  standard  notation 
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the  normal  vector  would  point  out  of  the  fluid;  In 
aeronautical  applications  It  Is  often  written  with  the 
normal  pointing  Into  the  fluid. 


(1.10) 


c/s^ 


-  /  // 

^JJP  p. 

-A 

Where  IP)  Is  the  outward  normal  to  S  into  the 
fluid  at  the  point 

coordinates  of  a  point  In  the 
fluid.  '7‘>S  coordinates  of  a  point  on  the 

closed  surface  S. 

The  surface  S  Is  taken  to  be  the  "upper"  cind  "lower" 
surfaces  of  a  zero  thickness  wing,  given  by  some  function 

.  In  linearized  theory  for  thin  wlng- 
llke  configurations,  (l.e.  for  <<  1  ) 

we  satisfy  the  upwash  boundary  condition  on  the  mean  position 
of  the  nonplanar  surface,  f  • 

We  rewrite  equation  (1.10)  as  an  Integral  of  the  Jump 


In  p  and  across  an  open  surface 

(i.i; 


(1.11) 
p: 
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Ap  -  P'-p"' 

S'*"  Is  the  "top"  of  the  surface 
S'  is  the  "bottom"  of  the  surface 

n  is  the  normal  to  the  surface  positive  "upward". 

With  an  extremely  nonplanar  surface  It  Is  often  quite 
arbitrary  which  Is  the  top  surface  and  which  is  the  bottom 
surface;  for  planar  wings  the  top  and  bo\;tom  are  more  easily 
defined. 

Pig.  (1.1)  shows  a  typical  nonplanar  lifting  surface 
configuration. 

The  nondlmenslonal  linearized  rac^ientum  equation  for  the 
direction  normal  to  S  Is 

for  S'*" 

+  ^ 

_  -h 

(1.12) 

for  S“ 

S  hi  ^  X  ^ 
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If  we  consider  a  zero-thickness  surface  In  fxilly  wetted 
flow  the  normal  velocity  Is  continuous  across  the  surface. 


Therefore,  for  fully  wetted  flow  the  function 
Is  zero  on  the  surface. 

Equation  (1.8)  then  becomes 


% 


^  (1-13) 


--L  ^  fj-) 

t  p./ 


The  singularity  “  ^  I  )  is  called  a  pressure 

doublet.  This  singularity  causes  a  unit  pressure  Jump 
across  the  surface  at  the  point  x  y  ^  , 

=■  ^  A7  Is  the  normal  at  that  point. 

The  position  of  a  nonplanar  surface  performing  unsteady 
motions  Is  given  by  some  function  of  3 , ‘y  j  'f  ^ 

For  small  perturbations,  the  upwash  velocity,  normal  to 
the  mean  position  of  the  surface  Is  formed  from 


^  Ot 


(1.15) 


1C 


•nie  velocity  at  a  point  x,  y,  z  In  the  direction 
n(y,  z  (y))  Is  given  by 

(1.16) 

The  velocity  potential  (p(x,  y,  z)  may  be  found  by  the 
operation  Indicated  In  equation  (1.5) 


(1.17) 


(1.18) 


where 


The  velocity  boundary  condition  (1.15)  Is  satisfied  In  the 
limit  as  the  point  x,  y,  z  approaches  the  surface 
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X,  y,  z  (y). 

The  kernel  function  for  this  problem  is 


J-  L 


-  Oo 


X-J 


(1.19) 


Referring  to  Pig.  1.1,  we  write  the  following  expressions 

N  N 

for  ^  and 


5h 


(1.20) 


(1.21) 
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Since 


aind 


^fi±J 


then 


(1.22) 


(1-23) 


i  a. 


r 

Gs.  ^(p  Gi  +  S/t! 

(1.24) 


) 


The  complete  form  for  this  kernel  function  Is  given  in 
Ref.  (8).  It  is  reproduced  in  Appendix  11.  For  a  planar 
lifting  surface,  the  kernel  function  becomes 


'^(  y- 


^  ck^-u-v) 

/^rx) 
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which  agrees  with  references  (6)  and  (j).  With  the 
kernel  function  defined  by  equation  (1.19)>  the  upwash 
velocity  normal  to  the  mean  position  of  the  nonplanar 
surface  becomes 


(1.26) 


/>,y,  ::(y; )  = 


This  integral  equation,  relating  the  complex  an^litude 
of  the  unsteady  load  on  the  surface  to  the  complex  amplitude 
of  the  upwash  velocity  normal  to  the  mean  surface,  is 
solved  numerically  by  a  technique  of  assumed 
modes.  This  method  is  discussed  in  section  (1.4),  and  a 
few  numerical  results  are  given. 

1.3  The  Kernel  Function  for  a  Nonplanar  Hydrofoil  at 

Infinite  Froude  Number 

For  a  hydrofoil  traveling  beneath  a  free  surface  at 
high  velocities  such  that  the  Froude  number  ^ 

can  be  taken  equal  to  infinity,  the  free  surface  becomes  a 
surface  of  zero  perturbation  pressure.  The  infinite  Froude 
number  approximation  neglects  the  effects  of  gravity  waves 
Induced  by  the  motion  of  the  foil.  In  Chapter  2,  these 
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effects  are  Included  and  numerical  results  show  that»  at 
least  for  foils  of  Infinite  aspect  ratio,  a  Proude  number  of 
10  Is  essentially  equivalent  to  Infinity  for  the 


prediction  of  loads  due  to  unsteady  motion. 

The  singularity  which  causes  a  unit  lift  at  the  point 


In  the  direction  as  well  as 

causing  no  perturbation  pressure  on  the  free  surface,  located 


at  z  -  0,  has  the  following  form: 


J.  A  /  .  '  . 

(1.27) 


where 


2  =  Cos 


(1.29) 


It  is  convenient  to  think  of  this  singularity  as  a 


pressure  doublet  with  axis  In  the  normal  direction  on  the 
foil  It  elf  and  a  pressure  doublet,  of  the  same  strength 
located  at  the  Image  point  with  axis  In  the  Image  normal 
direction.  Thus  a  foil  traveling  beneath  a  free  surface 
at  Infinite  Froude  number  Is  represented  by  a  foil  and 
Its  Image,  equally  loaded. 

The  kernel  function  for  this  problem  representing 
the  upwash  normal  to  the  foil  at  the  point  x,  y,  z,  due 
to  this  doublet  pair.  Is  given  by 


#7,  i-/)  J-  ^  [a 


>r-3 


T.tf) 


For  completeness,  the  final  expression  for  this  kernel 
function  Is  given  In  Appendix  11.  The  expression  for  the 
kernel  function  Is  singular  only  at  the  line  y  «  q,  z  •  'f 
on  the  foil;  the  singular  behavior  of  this  kernel  function 
Is  the  same  as  that  of  the  kernel  function  for  a 
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nonplanar  surface  In  an  Infinite  fluid.  The  numerical 
integration  can  be  handled  by  the  same  technique  as 
those  developed  by  Watkins  (Ref.  (7))  for  the  planar  wing 
in  an  infinite  fluid. 

1.4  The  Singularity  in  the  Kernel  Function 

The  kernel  functions  in  equations  (1.19)  #  (1.25)  and 
(1.30)  become  singular  on  the  line  y  =  z(y)  « 
on  the  surface  S.  This  is,  of  course,  a  familiar  situation 
and  the  solution  to  the  apparent  problem  is  to  find  the 
proper  definition  of  the  Improper  integral.  A  full 
treatment  of  Improper  Integrals  which  appear  in  two- 
dimensional  theory  was  given  by  Mangier,  Ref.  (11). 

Mangier  gives  definitions  for  improper  Integrals  which 
arise  when  the  limit  of  a  potential  function  is  taken. 

For  example,  the  Cauchy  principal  value  integral  is  the 
limit  of  the  following  potential  function. 


L 
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The  function  ^  is  a  solution 

V. 

of  Laplace's  equation.  Taking  the  t  derivative  of 
(1.31)  defines  the  Improper  Integral 


Ci'(t)  -  j f  J'T'fi-  2(-t-r)^  )] 

V  vi-ro !  J  I  j 

Mangier  shows  that  the  proper  definition  of  this 
Integral  Is 


h> 


In  order  to  apply  Mangier 's  results  which  are  valid 
for  functions  which  satisfy  the  two-dimensional  Laplace's 
equation,  we  must  show  that  as  we  approach  the  singularity, 
the  flow  becomes  a  locally  two-dimensional  potential  flow 
In  the  cross  flow  plane. 

We  Isolate  the  region  near  the  line  y  =  ,  z  = 


18 


by  a  strip  ^  of  width  2£  ,  along  the  chord.  In  this 

region  we  express  the  kernel  function  in  coordinates 
normal  and  tangential  to  the  surface  at  that  point. 

(See  fig.  1.2) 


Figure  1.2  The  Neighborhood  of  the  Singular 
Point . 

For  6  small  compared  to  the  radius  of  curvature 
of  the  surface  the  kernel  function  Is  approximately 


"hr 


ikp\  -tx-i)) 
d  A 


(1.3^) 


For  small  n  and  t,  equation  (1.34)  can  be  written 
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■  C  k(K-3) 

^rr 


IX-31A 


-  K(^-*.)  f-J—  - 

I  (i,vrV  / 

(See  Appendix  11). 

Other  slngiilarltles  are  present  but  are  of  lower  order. 
The  upwash  boundai^  condition  Is  expressed  as  the 
limit  of  the  normal  velocity  as  the  control  point  approaches 
the  surface. 


(1.36) 


In  the  narrow  region  CT  ,  the  pressure  Jump 

can  be  expressed  by  the  first  term  In  the  Taylor  series. 


^  Ap  (9, *  oce) 


(1-37) 
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The  form  of  equation  (I.36)  shows  that  as  we  approach 
the  line  y  *  ^  ,  z  «  ,  the  kernel  function  becomes 

a  two-dimensional  potential  function  in  the  cross  flow 
direction. 

The  integral  of  this  function  over  7^  and  the  limit 
as  n  0  has  been  shown  by  Mangier  to  be  properly  defined 
by  equation  (1.33)*  Thus  in  three-dimensional  potential 
flow,  the  singularity  of  the  kernel  function  for  the  upwash 
due  to  a  pressure  doublet  is  properly  handled  as  a  Mangier 
integral.  The  standard  notation  for  this  integral  is 


^(1.38) 

'  s-<r 


21 


1.5  Numerical  Method  for  a  Nonplanar  Fully  Wetted  Lifting 

Surface. 

The  Integral  equation  for  a  nonplanar  lifting  surface 
(1.38)  is  solved  by  essentially  the  method  develop  d  by 
Watkins  et  al.  Ref.  (7)»  for  the  planar  wing.  The  unknown 
load  distribution  p  on 

the  surface  is  represented  by  a  series  of  assumed  modes  with 
Unknown  coefficients.  These  coefficients  are  determined 
by  satisfying  the  upwash  velocity,  equation  (I.38),  at 
control  points  on  the  surface. 


AT  Cx,rt, 

' (1.39) 


where  y^^,  Zj^  are  coordinates  of  the  kth  point  and 
C  3 ith  assumed  mode. 

The  integrals  in  equation  (1.39)  are  evaluated  numerically 
by  Gaussian  quadrature  except  for  a  thin  strip  of  width  2  £ 
about  the  line  y  =  ,  z  =  ^  on  the  si.U’face.  If  the 

width  2  is  small  compared  to  the  radius  of  curvature 
of  the  surface  at  this  point,  the  method  developed  by 
Watkins  to  account  for  the  singularity  can  be  applied 
directly.  The  neighborhood  of  the  point  y  *  ,  z 

la  shown  in  figure  1.2.  We  define  a  coordinate  ^  to 
be  arc  length  measure  from  y  =  ^  ,  z  »  .  The  result  of  the 

intergratlons  in  this  strip  which  are  done  by  Gaussian  quadrature 
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Is  expanded  In  powers  of  S. 


(l-40) 

/ Api  ^  Cc,^  r** 

^  />» »  0 

then 


^  c/p' 

J.f.  7-=- 

If  the  surface  Is  synnnetrlc  in  y  and  If  the  upwash 
velocity  Is  purely  symmetric  or  purely  antisymmetric 
in  y,  the  function  L  =,  will  be 

symmetric  or  antisymmetric  in  The  boundary  condition 
in  this  case  need  be  satisfied  at  control  points  on 
only  half  of  the  surface. 

The  assumed  load  series  used  by  Watkins  for  a  plane 
wing  is 


H  M 


(1.4a) 


The  surface  is  located  such  that  the  root  semichord  is 
between  ^  *  -1,  the  leading  edge,  and  ^  **  +1>  the 
trailing  edge. 
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The  ^  functions  are  chosen  In  the  following  way. 

Let 

5  -  -  b(^)  ^  y-  ^4^-  (^J  ^ 

where  b(n)  Is  the  semichord  at  span and  3^^-^  Is  the 
leading  edge  coordinate  at  span 


=r  Cs€  f  /z) 

-  fi^  eci)) 


(1.44) 


The/ (3)  functions  satisfy  the  Kutta  condition  at 
the  trailing  edge  and  the  /{$)  function  has  a  square  root 
singularity  at  the  leading  edge.  Analytical  results  for 
a  flat  plate  In  two-dimensional  flow,  and  for  a  wing  of 
circular  plan  form  In  three  deraenslons,  show  a  square 
root  singularity  at  the  leading  edge. 

The  functions  for  symmetric  upwash  are 


The 


/xh)  //-/V/s;^ 

^^(/y)  functions  for  antisymmetric  upwash  are 


{1-45) 


where  is  the  semi  span.  The  /«(7)  functions  satisfy 


(1.46) 


the  condition  that  the  lift  drops  to  zero  at  the  wing 
tips.  They  also  are  the  functions  used  In  lifting  line 
theory.  In  general,  one  should  pick  a  set  of  functions 
to  have  properties  as  close  as  possible  to  the  actual 
lift  distribution.  In  this  case,  control  point  location 
Is  less  critical  and  convergence  Is  possible  with  fewer 
modes . 

For  a  nonplanar  surface,  we  could  use  this  set  of 
functions  or  replace  the  spanwlse  functions  with  a  set 
more  appropriate  to  a  partlcvJar  configuration.  For 
surface  piercing  foils.  It  is  not  likely  that  the  lift 
drops  to  zero  at  the  water  surface  with  Infinite  slope. 

For  Intersecting  surfaces  the  load  does  not  drop  to  zero 
at  the  comer  and  there  Is  the  additional  constraint  that 
the  sum  of  the  pressure  Jumps  at  the  Intersection  Is  zero. 

The  result  of  these  techniques  give  a  set  of  complex 
linear  algebraic  equations  for  the  unknown  coefficients 


of  the  L. p 


modes . 
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Once  the  b;s  have  been  found,  the  pressiire  Jisnp 
distribution  on  the  surface  is  known.  Prom  this,  force 
and  moment  coefficients  for  the  configuration  may  be 
calculated. 

The  numerical  technique  has  been  programmed  for  the 
IMB  7090.  In  practice,  for  simple  configurations,  calculations 
with  9,  16,  and  25  control  points  have  given  good  results. 

Time  required  for  a  calculation  is  under  five  minutes. 

As  an  example  of  the  method,  calculations  were  done 
to  compare  with  the  experiments  of  Abramson  and  Ransleben, 

Ref.  (12).  For  these,  the  unsteady  lift  and  moment 
coefficients  at  several  chordwise  sections  on  the  span 
of  an  aspect  ratio  5  foil  in  bending  and  torsion 
oscillations  were  measured.  Pig.  I.3  to  1.7  show  numerical 
results  for  unsteady  lift  coefficient,  magnitude  and 
phase,  along  the  span  due  to  bending  at  reduced  frequencies 
K  «  .6,  .8,  1.2,  2.0.  These  calculations  were  carried 
out  with  S,  16  and  25  control  points.  Results  are  in  fair 
agreement  at  low  reduced  frequency;  agreement  is  poor  at 
higher  reduced  frequencies . 

This  is  perhaps  due  to  a  failure  of  the  Kutta 
condition  at  high  reduced  frequencies.  This  is  discussed 
by  Ashley,  Wldnall  and  Landahl  in  Ref.  (9)  and  some 
empirical  modifications  to  the  numerical  method  are  suggested. 
Calculations  were  also  performed  for  a  planar  wing  below 
the  free  surface.  These  results  were  presented  in  Ref.  (8), 


THEORY  (16  PTS.) 
EXPERIMENT(K=2  0) 


FIG.  1.6  COMPARISON  WITH  EXPERIMENT  -  SPANWI SE 

DISTRIBUTION  OF  LIFT  DUE  TO  BENDING  (  k  =  2.0) 


30 


CHAPTER  II 

THE  EFFECT  OF  GRAYITY  WAVES  ON  THE  UNSTEADY  LIFT  AND 
MOMENT  ON  AN  OSCILLATING  HYDROFOIL 

2.1  Introduction 

An  intriguing  question  in  the  study  of  unsteady  hydro¬ 
foil  motions  is  whether  the  gravity  waves  Induced  by  the  motion 
Influence  the  loading  to  any  great  extent.  To  answer  this  ques¬ 
tion  we  formulate  the  ’’exact  linearized"  Integral  equation  rela¬ 
ting  the  loading  on  an  oscillating  foil  near  the  free  surface 
to  the  upwash  on  the  foil.  In  this  problem,  we  represent  the 
loading  by  a  distribution  of  singularities  on  the  foil  surface 
which,  in  addlt.lon  to  producing  a  unit  lift  at  a  point,  have  a 
perturbation  pressure  field  which  satisfies  the  free  surface 
boundary  condition.  The  Integral  equation  will  then  be  solved 
by  a  method  of  assiuned  modes  similar  to  that  described  in  Chap¬ 
ter  I . 

V'e  locate  our  coordinate  system  so  that  the  mean  posi¬ 
tion  of  the  free  surface  is  the  x,y  plane  at  z  =  0.  The  free 
stream  velocity  is  in  the  positive  x  direction.  As  seen  in 
this  coordinate  system,  the  free  surface  will  take  some  shape 
P  =<^(x,y,t).  The  boundary  condition  on  the  free  surface  is 
that  the  perturbation  pressure  is  equal  to  zero.  (See  ref.(f^/)). 

The  linearized  Bernoulli  equation  for  an  incompressible 
Irrotatlonal  flow  in  a  gravity  field  is 
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lA  t  p  2^.  -o 

<st 


(2.1) 


where  the  variables  have  been  nondlmenslonallzed  with  respect 
to  free  stream  velocity^,  density p,  and  root  semichord 

P  Is  the  Proude  number  (P  » 

^  Is  the  perturbation  velocity  potential 

A 

P  Is  the  perturbation  pressure 
Applied  at  the  free  surface  ^  (x,y,t)  this  equation  becomes 

'^t  ^x. 

The  vertical  velocity  of  a  fluid  particle  on  the  free  surface 
is  related  to  the  velocity  potentla^^.  by 


(2.3) 


A  combination  of  equation  (^.2)  and  (2.3)  gives  the 

I 

free-surface  boundary  condition 


{2.U) 


In  the  linearized  problem,  this  boundary  condition  Is 
applied  on  the  mean  position  of  the  free  surface,  z  »  0. 

The  same  boundary  condition  must  also  be  satisfied  by 
the  perturbation  pressure. 


rz 

For  simple  harmonic  motion  with  reduced  frequency  k 


(2.5) 


(2.6) 


The  complex  amplitude  of  F  satisfies 


(2.7) 


If  the  Froude  number*  is  infinite,  this  boundary  con¬ 
dition  reduces  to  zero  pressure  perturbation  on  the  free  surface. 
As  was  discussed  in  Chapter  I,  the  singularity  which  causes  a 
unit  lift  at  a  point  and  causes  no  perturbation  pressure  at 
Z  =  0,  is  a  pressure  doublet  plus  its  image. 

To  find  the  singularity  which  causes  a  unit  lift  at  a 
point  on  the  foil  and  in  addition  produces  a  pressure  field  which 
satisfies  the  full  free  surface  boundary  condition  is  more  com¬ 
plicated.  For  two-dimensional  flow,  analytical  expressions  for 
the  kernel  function  can  be  found  using  Fourier  transform  tech¬ 
niques.  The  problem  in  three  dimensions  can  be  reduced  to  a 
single  Integral  which  could  be  done  .numerically. 

When  the  upwash  on  the  foil  due  to  these  singularities 
has  been  found,  the  problem  may  be  expressed  as 


where  the  integration  is  over  the  surface  of  the  foil.  All 
effects  due  to  the  free  surface  are  contained  in  the  kernel 
function. 


2.2  The  Kernel  Function  for  Two  Dimensional  Flow  at  Arbitrary 
Froude  Number 

The  complex  amplitude  of  the  perturbation  pressure 
fields  for  a  lifting  element  at  z  -  -d  satisfies  Laplaces 
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equation  (2.9)  with  boundary  conditions  (2.10)  -  (2.13) 


7Y-0  (2.9) 

« 

=  °  ^  <2.10) 

pc^;d')  - p(y^.J*)  =  /cx-O  (a. 11) 

^  -J*)  =(5  (a.ia) 

P(3L,-o^)  -O  (2.13) 

The  complex  amplitude  of  the  upwash  velocity  due  to 

this  pressure  field  Is  given  by 


f  /'*  CK(x-Cx-i)) 

In  addition  to  these  conditions,  we  must  impose  a 
radiation  condition  requiring  that  all  waves  we  consider 

j 

have  been  generated  by  the  motion  of  the  foil.  This  may  be 
done  by  letting  the  reduced  frequency  k  have  an  infinltessimal 
negative  imaginary  part  so  that  at  time  t  ■  -o*,  there  was  no 
motion. 


The  derivation  of  the  kernel  function  is  similar  to 
that  of  ref.  (8),  We  apply  Fourier  transfonns  in  x  according 
to  the  definitions 


f>  (:d  ^ -2  ) 


e 


-«-'5 


y 


(2.15) 

(2.16) 


3^ 


/%• 

(2.17) 

oc 

y^o 

Equations  (2.8)  to  (2.13),  In  the  transform  variables  and 

become 

-s' £  =o 

(2.19) 

-a*st£ * -y.  ^ -o 

2-0  (2,20) 

Z' (^,-d)  -  rZ (^,-d)^  e~‘^^ 

(2.21) 

<i_2~  (i,-J)  -  <^Z^(',-J)  r' 

jy  Jz 

(2.22) 

Z  (s,  -<=°)  =  o 

(2.23) 

ii  ^  J-  dj 

dz,  (k*s) 

(2.24) 

where  J?  Is  the  solution  for  and  ^  Is  the 

solution  for 

z  -T^-  d  , 

•me  solutions  to  equation  (2.12)  are 

^  isiz.  -mz 

=  /i  e  -  ^  e 

(2.25) 

/s/  z.  -  “IV  Z 

Z  ■-  h  e.  B  e 

(2.26) 
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Condition  (2.23)  requires  *=  0.  The  functions  ABA 
are  found  from  equation  (2.20  -  2.22).  The  resulting  expressions 
for  and  •C  are 


•(2^  /s/  -  L  SJ 

e 

'o/j/s/  -t  sj  V 


(2.27) 

(x-tj)fsi  -tSjy 

h  j—  <2.  el  ^  ) 

j*-.;  /s^ 

(2.28) 

The  kernel  function  can  be  found  from  either  (2.27)  or  (2.28) 
by  using  equation  (2.24)  and  (2.1') 


--  !  fTir  "'^/(a.30) 

The  first  term  In  equations  (2.27)  and  (2.28)  gives  the  kernel 
function  for  a  pressure  doublet  In  an  Infinite  fluid.  Thus 

00  —6/s/  ^  CsC^c -"i) 

^  i/js^  e  c^s  j 

We  use  the  fact  that  k  has  a  small  negative  Imaginary 
part  to  move  the  pole  off  the  axle  of  Integration.  Hence 


;here  results 


(2.32) 


i<  O' ‘m)  )7 


>6 


where 


(2.33) 


Ec(y^)  = 

^'oo 

is  the  complex  exponential  integral  with  branch  line  along  the 
positive  real  axis. 

The  prop''.'  definition  of  an  Integral  over  the  general¬ 


ized  function 


8  the  Cauchy  principal  value. 


A-  (X'3>)  *6^ 


'  ?  =  I  di 
j  /  r>-3) 

6i» 


(2.34) 


The  first  term  In  equation  (2,32)  contains  the  strongest  singu¬ 
larity  In  the  kernel  function.  The  second  term  has  only  a 
logarithmic  singularity  and  the  additional  terms  due  to  free 
surface  effects  are  not  singular. 

For  Infinite  Proude  number,  the  kernel  function  Is  Just 


e 

f  f 


l/i/ 


-6lsU^sCK^t)^ 

£  (ds  /  /  ! } 


-  09 


-XJh/-h  LiCx-S) 

e  </s 

{2-35) ^ 


which  gives 


s:  —  -Z- 
2jr 


2}r 


e  I 


(2.36) 


wnere 
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Zd  4- L  /  V  -  3  9 
c  /x  -?j) 


(2.37) 

(2.38) 


It  Is  Interesting  to  note  the  asymptotic  form  of 
the  kernel  function  for  large  values  of  x  -  3.  The  asyrr^* 
totlc  representation  of  the  complex  exponential  Integral  Is 


E.  ~  ( ^)  ^ 

The  lowest  ox\ier  terras  In  equation  (2,3^)  are  thus 


(2.39) 


A( 


k  f  <  IL  ( 


zrr  L  ^  I 

viJ 


(  v-'S) 
4 


)] 


-2  -  c‘k 


{2, ho) 


■S  rr  ~  (  k’-t) 

The  behavior  far  downstream  from  the  disturbance  is 

that  of  a  wave  moving  In  the  positive  x  direction  with  unit 
velocity. 


/<.  '<  (  I  t 


( X-f)  \ 


e' 


k(y~ty 


/x-f/J  ^  (2.“!) 

No  wave  exists  upstream,  'Hils  represents,  of  course,  the  shed 
wake  vortlclty,  which  Is  carried  downstream  by  the  free  stream. 
Later  on,  we  shall  use  these  same  Ideas  to  find  the  gravity 
waves  produced  by  the  motion.  At  this  point,  we  will  note  that 
the  poles  Inflnltesslmally  close  to  the  axis  of  Integration  in 
eq.  (2,33)  give  rise  to  waves  which  travel  tjpstream  or  down¬ 
stream,  If  the  poles  had  been  located  at  a  finite  distance 
from  the  real  axis,  the  waves  would  be  damped. 
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4 


/v-^/  '  (2.^2) 

The  complete  expression  for  the  kernel  for  finite  Froude 
number  Is 


<S- 


fx-i. '  ^ 


p  r  f  isCn-rj  --'5' 

■=  ~  ^  L  i 


oo 


isC^t  -i)-  2dlSi 
^  _1  I L  e.  1^1 


9C 


^  -9  0  I  ^iT i  ^  (TTk) 

/  ^  1^1  )  ds: 

(  ISI  /  (2.' 


(^t):)  ^  -  ISI  /  (2.43) 

•» 

where  K  Is  the  kernel  function  for  a  pressure  doublet  In  an 
infinite  fluid  and  is  the  correction  term  due  to  free  surface 
effects.  This  correction  term  may  be  expressed  as 

^  -frr  I  J  (ifi)  '  FVci-n)‘*s  / 

'  /  ,  v-a  S  2sd  -  t  V  } 

^  (  SJldlsS)  e  ?  ,^\n) 

The  first  Integral  has  poles  at 

-  A 

Si  ^  / 

f  z 

where  f=I^F  .  'Riese  are  all  located  infinitesimally  close 
to  the  real  axis  along  the  path  of  integration, 

The  second  Integral  has  poles  at 

=■  k 


(£'.‘15) 

/  V-  *  2F  ^ 

(£'.‘i6) 

■‘rj  -^7^r/wJ 

(£'.■'7) 

k( I- 


a/  -2f 


h-4^) 


(2.49) 
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5^*  k( /'  ff  -  j/  j  (8. 30) 

The  location  of  these  poles  depends  on  the  value  of  f.  Thus, 
for  they  are  located  Infinitesimally  close  to  the 

real  positive  axis.  For  they  are  complex  conjugates 

in  the  second  and  third  quadrants.  In  the  case  of  ,  they 
are  complex  conj\igates  in  the  first  and  fourth  quadrants.  These 
poles,  and  their  significance  for  the  gravity  waves  generated 
by  the  motion,  have  been  discussed  by  Tan,  Ref,  (12),  and 
Kaplan,  Ref.  (13)  and  (l4). 

A  root  locus  diagz*am  is  shown  in  fig.  2.1  for  constant 
F,  and  increasing  k, 
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2.1 


Location  of  the  poles  of  equation  (2,44)  for 
constant  P,  increasing  k.  (f  «  kP^) 
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The  Integrands  In  equation  (2,44)  may  be  separated  Into 
partial  fractions 
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where 


didfl  (i+  f')i 


(2.52) 

(2.53) 
(2.5'i) 
(2.55) 


The  first  Integral  in  equation  (2.51)  l8  the  kernel  function 
for  a  simple  steady  pressure  doublet  at  the  image  point.  The 
remaining  Integrals  can  be  expressed  as  complex  exponential 
Integrals  with  or  without  a  correction  term  depending  on  whether 
the  path  of  integration  of  equation  (2,5l)»  after  transfbrmatlon 
into  the  exponential  Integral  form, equation  (2,33)#  crosses 
the  branch  line  of  this  function, 

'Oils  depends  on  the  position  of  the  poles  in  the  complex 
plane  which  is  a  function  of  f.  (f  *=  kP^) 


For  the  kernel  function  Is 
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(2.60) 
(2.61) 


For  large  values  of  (y  ”  5)»  asymptotic  form  of  the  kernel 


for  f  is 
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Since  for  ^  ^  are  real  and  negative,  we 

can  draw  the  conclusion  that  the  disturbance  produces  four 
gravity  waves,  in  addition  to  the  wake.  The  s^^  and  s^  waves 
travel  downstream  behind  the  disturbance,  the  s^  wave  travels 
upstream  ahead  of  the  disturbance,  and  the  S|^  wave  travels 
upstream  behind  the  disturbance. 

In  the  limit  as  the  s^  and  s^  waves  have  a 

group  velocity  idilch  Is  equal  to  the  velocity  of  the  foil 
relative  to  the  stationary  fluid,  (see  Ref,  14)  This  Is  a 
singular  point  for  the  linearized  problem,  in  some  ways  analo 
gous  to  the  Mach  number  unity  problem  for  linearized  compres¬ 
sible  flow. 

For  f  '/f  ,  the  kernel  function  differs  only  In  the 
contribution  from  the  poles.  Thus  for  Kj  V': 
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whereas  for  f  >  ?  the  kernel  function  Is 


y-i4,  £i(<ij  i- ii4  e'^" 

^rr  ^  ^ 

^  ^  ElLcl^) 

(2.64) 

The  asymptotic  form  of  the  kernel  function  for  J  ^  f  shows  that 
the  gravity  waves  from  the  s^  and  s^  poles  are  damped  far  from 
the  origin.  Only  the  gravity  waves  due  to  the  s^  auid  s^  poles, 
plus  the  unsteady  wake,  extend  downstream  from  the  disturbance. 


2.3  Numerical  Solution  and  Results  for  Arbitrary  Proude  Number 

The  Integral  equation  relating  the  complex  auaplltude 
of  the  upwash  and  the  load  distribution  on  a  two-dimensional 
oscillating  hydrofoil,  traveling  beneath  a  free  surface  Is 

-f  / 

Where  K(x-^)  «  K*  +  as  derived  In  section  2.2.  As  In  the 


three-dimensional  case,  the  kernel  function  is  the  upwash  at 
a  point  X  due  to  an  oscillating  pressure  doublet  at  a  point 
nils  function  contains  the  effects  of  the  unsteady  wake,  the 
surface  waves  generated  by  the  motion,  and  the  depth  of  the 
foil  below  the  free  surface. 

A  numerical  solution  of  this  equation  has  teen  pro¬ 
grammed  for  the  IBM  7090.  The  method  which  was  used  Involves 
a  series  of  assumed  modes,  in  a  manner  similar  to  that  described 
in  Chapter  I,  One  assumes  a  series  of  pressure  modes,  each 
of  which  satisfies  the  Kutta  condition,  and  solves  for  the  co¬ 
efficients  of  these  modes  from  the  known  upwash  distribution. 
Once  the  coefficients  of  the  pressure  modes  are  known,  the 
unsteady  lift  and  moment  can  easily  be  calculated. 

The  foil  is  located  on  the  J  axis  from  3  =  +1  to  S  = 
the  semlchord  thus  being  adopted  as  a  reference  length.  The 
assumed  lift  functions  are  chosen  as  follows; 
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3  =  -  cos  ©  (2.68) 

The  are  the  unknown  complex  coefficients.  With  this  assump¬ 


tion  for  the  lift  distribution,  equation  (1.1)  becomes 

•7?^  4. 
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(2.69) 


Define  a  function 
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In  terms  of  the  , 
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the  Integral  equation  becomes  a  com¬ 


plex  matrix  equation 


(2.71) 


where  Is  the  known  complex  upwash  evaluated  at  a  preassigned 
set  of  control  points  x^.  It  would  be  possible  to  solve  for 
N  by^’s  by  fitting  the  upwash  at  N  points.  Then  would 

be  a  square  matrix.  Instead,  the  upwash  Is  "fit”  by  a  mini¬ 
mum  meaui-square  error  criterion  at  more  points  than  the  number 
of  assumed  modes.  This  device  makes  control-point  location 
less  critical. 

Results  of  an  analyslr  to  choose  the  complex  coefficients 
for  a  "best"  fit  have  led  to  the  following  procedure;  One 
writes  a  linear  set  of  equations 


MxW 


(2.72) 


where  M  <  N. 

Then  premultiply  both  sides  by  the  complex-conjugate  transpose 

0^“  fcch]  • 


(2.73) 


The  result  Is  an  Mth  order  set  of  linear  equations  for  the  M 
unknown  complex  coefficients.  This  procedure  minimizes  the  dot 
product  of  the  error  and  the  conjugate  of  the  error.  This 
method  Is  used  In  the  program  to  fit  the  coo^lex  vqpwash  v(x)  by 
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the  sum  of  complex  upwashes  due  to  the  assumed  lift  modes. 

In  the  present  program,  nineteen  downwash  points  are 
chosen  and  a  "best  fit"  is  made  with  four  or  five  pressure 
modes.  When  the  conq)lex  coefficients  have  been  found,  unsteady 
lift  and  moment  coefficients  are  computed  tyttie  familiar  formulae 


Cl- ^ir  ( k, /s.) 

rrU(  h.-b,)  (2-T5) 

The  program  for  infinite  depth  reproduced  the  Theodorsen 
results.  Ref,  (1),  to  four  significant  figures. 

Results  for  and  for  a  two-dimensional  hydro¬ 
foil  oscillating  in  heave  and  pitch  beneath  a  free  surface 
are  shown  in  Pig,  2,2. to  2,5. 

Definitions  of  lift  and  moment  coefficients  are 


(2,79) 


By  these  definitions «  has  a  quasi-steady  value 
oT  2'ir  at  infinite  depth. 

Results  are  shown  for  lift  and  moment  coefficients 
vs.  reduced  frequency  for  various  Proude  numbers  and  depth 
below  the  free  surface. 

These  results  indicate  that  the  coefficients  predicted 
for  Proude  numbers  of  about  10  are  equivalent  to  those  predicted 
for  infinite  Proude  numbers.  Thus  for  hydrofoil  operation  in 
typical  ranges  of  P,  the  effects  of  free  surface  waves  are  not 
significant. 

The  behavior  of  these  coefficients  near  the  singular 
point  kP^  =  1/4  is  Interesting  from  a  mathematical  point  of 
view  but  seems  to  be  of  limited  practical  Importance  in  prob¬ 
lems  of  hydrofoil  flutter. 

2.4  The  Kernel  Punction  in  Three  Dimensions 

To  determine  the  form  of  the  kernel  function,  we  make 
\i8e  of  Fourier  transforms  in  the  x  and  y  direction.  Definitions 
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of  these  transforms  are 


,5^ /y  ^  ’Js 
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■“’  (2.86) 


The  relation  between  the  and  ^  transforms  is 

^  (2.87) 

Laplace's  equation  for  (p  aind  p  becomes  an  ordinary  second- 
order  differential  equation  for  f  and  ^  ,  with  i  as  the 
Independent  variable.  The  transform  variables  appear  as 
parameters 


dll  -  fs  dF 

ct  Z.^ 

dli.  _  ~  ^ 
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The  free  surface  boundary  condition  for  Is 


(2.88) 


(2.89) 


-  Z  ^  y-2.  ^  ■=  <^  <S’^- 


(2.90) 


Also,  as 
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(2.91) 


The  basic  equations,  as  Just  derived,  hold  for  general 
surfaces  at  arbitrary  orientations  below  the  free  sui’face.  To 
achieve  clarity  in  presentation,  we  shall  now  restrict  our¬ 
selves  to  lifting  surfaces  oriented  parallel  to  the  free  surface. 

We  located  the  planar  lifting  surface  at  z  «  -d  below 
the  free  surface.  The  kernel  function  for  this  surface  has  the 
property  that  across  z  ^  -d,  P(x,y,2)  experiences  a  Jump  of 
unit  strength. 


d  )  H-y) 


(2.92) 


The  Jump  conditions  on  the  transform  of  P(x,y,z)  are 


'  2.r 

y.  (2-95) 

where  ^  is  the  Fourier  transform  for  z<  -d  and  ^  is  the 
Fourier  transform  for  z  >  -d. 

For  a  zero  thickness  lifting  suiface,  the  upwash  is 
continuous  across  z  -  -d.  This  condition  becomes  a  condition  on 
of 


O'Z' 


(2.94) 
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The  two  solutions  to  equation  (2*80)  are 


Me 


2'  -  C 


z 


(2.95) 

(2.96) 


To  solve  for  A,B  and  C  we  apply  the  boundary  conditions 
(2.90)  at  z  «  0,  and  the  Jump  conditions  (2.93)  (2.94)  at 

z  -  -d.  In  order  to  insure  that  O  ,  as  z  -r  -  c»  ,  the 
real  part  of  the  radical  In  X  must  be  taken  as  positive.  Appli¬ 
cation  of  these  conditions  gives 
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For  an  oscillating  pressure  Jump  at  z  «  ~d  in  an  Infinite  fluid 
equation  (2.97)  reduces  to: 
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(2.99) 
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For  an  oscillating  pressure  Jian^)  at  z  =  -d  In  a  fluid  at 

.ii)  \riO 

infinite  Froude  number,  we  have  (2 .99)*  (2^01)  plus  a  correc¬ 
tion  term  due  to  a  positive  image  located  at  z  =  +d. 


^'r  *  ^ 
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The  normal  i^pwash  on  the  planar  surface  is 

^Cx,^)  -  (x ^4^,  d  i  k) 

az. 


(2.103) 


The  expression  for  the  velocity  potential  transform  is 


(2.10't) 


The  upwash  due  to  a  unit  oscillating  lift  located  at  x  -  S, 
y  ,  z  =  -d  on  a  planar  surface  parallel  to  the  free  sur¬ 
face  is 

0-0  ^  (S  xf  f  /S 

i^-d  [  ilr  [f  ^  (2.105) 

The  resulting  kernel  function  valid  for  z  -d  is  found  by 

applying  equation  (2.105)  to  the  form  of  in  equation  (2.97) 


O0 
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(isCx-s) 

(2.106) 


The  first  integral  in  equation  (2.106)  is  a  well  known  expres¬ 
sion  for  the  upwash  due  to  an  oscillating  pressure  doublet 
(see  Chapter  I). 

We  will  discuss  the  second  integral  and  propose  a  method 
for  numerical  evaluation  and  an  interpretation  of  the  results. 

The  complete  correction  Including  the  image  is 

<s*o 

^  iff’-  jf  ■sx/o  • 

- 

/  ]  ds 

'  /  \l=^(Stk)'^-  t's^F/S^/  '  (2-107) 


A  change  of  variables  which  gives  a  striking  comparison  between 
the  two  and  three-dimensional  kernel  functions  is  s  «  s,  3  »  ts 
where  t  is  a  real  number.  Properties  of  this  transformation 
are  worked  out  in  A-ppendix  ill.  General  results  are  that  if 
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With  this  transformation,  the  expression  for  the  kernel  due  to 
a  pressure  jump  In  an  Infinite  fluid  Is 


/;v 
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e'xpC  ^ dsJv I (2.110) 


This  may  also  be  written 
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The  e2q)resslon  for  the  equivalent  two -dimension -1  kernel 
function  for  a  pressure  jump  In  an  Infinite  fluid  is 
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Comparing  "V®®  and  ')/ 


we  notice  that 
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We  may  write  a  similar  expression  for  the  upwash  contribution 
of  the  free  surface  Including  the  Image.  Equation  (2.107)  becomes 

o»o  o* 

•  exp(  csu-i^  F7r)di 
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Ab  discussed  in  section  2,2,  the  equivalent  expression  for  the 
two-dimensional  kernel  is 

CJil  e  /  \ 

^ )  (2.115) 

The  expression  for  the  three-dimensional  kernel  becomes 
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An  analytic  closed  fom  expression  for  was  de¬ 
rived  in  section  2.2.  The  advantage  of  an  expression  like 
(2.116)  is  that  radiation  conditions  carefully  worked  out  for 
the  two-dimensional  kernel  can  be  applied  to  this  form  of  the 
three-dimensional  kernel.  This  gives  a  unique  definition  to 
the  various  contour  Integrals  vdilch  are  evaluated  in  solving 
for  this  three-dimensional  kernel  function.  Equation  (2.116) 
in  addition  to  settling  questions  of  incoming  and  outgoing 
waves,  also  gives  a  numerical  technique  for  evaluating  the 
three  dimensional  kernel  function.  Convergence  properties 
can  be  checked  by  examining  the  asymptotic  e^ansions  of 
^  ^20  for  large  T, 


The  final  result  for  the  three- iimensional  kernel 
function  is  reminiscent  of  an  ejqsansion  in  plane  waves.  First 
we  note  that  the  "apparent  depth"  of  the  equivalent  two-di¬ 
mensional  disturbance  grows  as 


The  apparent  Proude  n\unber  of  the  Individual  two-dimensional 
contributions  decreases  as 
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(2.118) 


For  a  given  value  of  ITl  ,  waves  come  from  both  the 
redirections.  A  better  form  to  demonstrate  the 
equivalence  of  the  and  directions  would  be 
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For  hydrofoils  operating  at  practical  ranges  of  Froude  numbers, 
it  16  felt  that  the  ntimerlcal  evaluation  of  these  Integrals 
would  not  be  practical.  In  the  two-dimensional  finite  Froude 
number  solution,  results  of  F  =  10  were  Indistinguishable  from 
infinite  Froude  number  results  for  all  practical  purposes. 

The  computation  time  required  to  carry  out  the  evaluations  of 
the  kernel  function  necessary  for  a  determination  of  the  loads 
of  a  three-dimensional  hydrofoil  at  flr»lte  Froude  number  could 
not  be  Justified. 

Extension  to  non-planar  lifting  surfaces  can  be  made 
by  noting  that  the  nature  of  the  singularity  which  gives  a 
pressure  Jump  across  a  surface  is  the  normal  derivative  of 
a  pressure  source,  the  pressure  doublet.  Taking  advantage 
of  this,  we  derive  the  expression  for  a  pressure  source  below 
a  free  surface  and  the  expression  for  the  upwash  it  produces. 


We  then  take  a  derivative  In  an/  direction  we  choose  to  find 
the  kernel  function  for  an  arbitrary  nonplanar  lifting  surface 
The  problem  statement  is  essentially  that  for  the 
planar  kernel  function.  The  jump  conditions  will  be  those 
associated  with  a  pressure  source  which  when  differentiated  in 
the  direction  gives  the  unit  pressure  jump  for  a  planar 
foil  parallel  to  the  free  surface. 

The  jump  conditions  on  pressure  across  this  pressure 
source  located  atz-J^,x«^  ,  y  m  are 

a  'I,  f)  =  ^^XKLxzx) 

In  the  transform  ven^lable  ^  the  jump  conditions  are 
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(2.125) 


These  equations  together  with  equations  (2.88)  to  (2.91)  give 
the  solution  for  P"*"  and  P"  for  this  pressure  source  located 
at  z  .f 

-cxp(-[^s-  ‘■M  -  (z-f))  ! 


(2.124) 
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To  find  the  normal  upwash  at  a  point  x^y^z  of  a  surface 
due  to  a  unit  pressure  jump  at  a  point  of  a  surface, 

we  follow  a  procediire  similar  to  that  of  Chapter  I  for  non- 
planar  wings  in  an  infinite  fluid  or  at  infinite  Proude  number 
below  a  free  surface. 

We  derive  the  kernel  function  for  three-dimensional 
finite  Proude  number  by  the  following  operation  on  the  velocity 
potential  due  to  the  pressure  source. 


(2.128) 


where  as  In  Chapter  I 


Cos '/'(if) 

(2.129) 

-  -f  j/n  A 

(2.130) 

C^J, 

The  configuration  a  nonplanar  surface  is  sketched  4n  Fig.  1,1 

If  the  kernel  function  described  in  equation  (2,123) 
could  be  evaluated  either  analytically  or  numerically,  it  would 
replace  the  kernel  function  used  in  Chapter  I  if  free  surface 
effects  were  desired.  When  the  Proude  number  P  is  infinity, 
the  kernel  function  described  by  equation  (2,128)  reduces  to  the 
simple  pressure  doublet  plus  image  doublet. 


FIG. 2.4  MOMENT  COEFFICENT  (  I/4C)  DUE  TO  HEAVE  VS  k 
FOR  A  TWO  DIMENSIONAL  HYDROFOIL 
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CHAPTER  III 

UNSTEADY  LOADS  ON  SUPERCAVITATINQ 
THREE  DIMENSIONAL  HYDROFOILS 
3.1  Introduction 

A  hydrofoil  traveling  at  high  speed  vrlll  experience 
cavitation  If  the  pressure  on  the  foil  falls  below  some 
critical  value.  If  the  cavity  extends  over  the  complete 
upper  surface  of  the  foil,  the  foil  Is  said  to  be 
supercavltating.  The  pressure  in  the  cavity  may  be  the 
vapor  pressure  or  In  the  case  of  operation  close  to  the 
free  surface,  it  may  be  atmospheric  pressure.  The  latter 
condition  is  often  called  ventilation.  If  a  cavitation 
bubble  forms,  it  will  have  a  constant  pressure  on  its 
surface  equal  to  this  critical  value.  To  solve  math¬ 
ematical  problems  of  cavity  flows,  it  is  usually  assumed 
that  this  critical  pressure  is  given. 

Problems  of  finite  cavltating  bodies  such  as  vertical 
flat  plates  in  two-dimensional  steady  flow  have  been  solved 
using  complex  variable  techniques.  A  review  of  these 
problems  and  the  mathematical  techniques  used  Is  given  by 
Gilbarg,  Ref.  (26).  Some  new  developments  in  two 
dimensional  cavitation  theory  are  contained  In  papers 
of  Ref.  (27). 

For  thin,  two-dimensional  hydrofoils  In  steady  and 
unsteady  motion,  linearized  theories  have  been  studied 
by  Tulin  (28), Woods  (29),  Parkin  (30),  Tlmman  (31) 
and  Guerst  (32). 
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In  his  doctoral  desertatlon,  Guerst  has  discussed 
the  relation  between  the  two-dimensional  linearized  theory 
for  a  cavltatlng  wedge  and  the  nonlinear  models  for  a 
cavltatlng  wedge,  the  Rlabouchlnsky  model  and  the  re¬ 
entrant  Jet  model.  He  shows  that  "linearized  cavity 
theory  Is  a  first  order  approximation  to  both  nonlinear 
models  for  small  wedge  angle." 

The  two-dimensional  linearized  theory  for  steady 
flow  have  shown  good  agreement  with  experiment  for  thin 
supercavltatlng  hydrofoils.  Therefore,  we  feel  that 
a  three-dimensional  linearized  theory  will  give  accurate 
predictions  for  the  loads  on  steady  and  unsteady  super¬ 
cavltatlng  foils  of  finite  span. 

3.2  Linearized  Theory  for  Three-dimensional  Supercavltatlng 

Hydrofoils 

The  coupled  linearized  Integral  equations  for 
cavltatlng  hydrofoils  can  be  put  Into  the  general  frame¬ 
work  of  lifting  surface  theory  ujlng  the  pressure  or 
acceleration  potential  and  Green’s  theorem.  For  linearized 
Incompressible  Irrotatlonal  flow,  the  nondlmenslonal 
perturbation  pressure  and  the  perturbation  velocity 
potentials  satisfy  Laplaces  equation. 
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(3.1) 

(3.2) 

In  the  application  of  Green’s  theorem  to  the  present 
problem  we  will  use  as  the  surface  S  enclosing  the  fluid, 
the  wetted  surface  of  the  foil  and  the  cavity  surface 
plus  a  surface  at  Infinity.  Defining  the  unit  normal 
n  as  pointing  Into  the  fluid  from  the  closed  surface  we 
then  obtain 


where 


(3.3) 


R  =  radius  vector  from  field  point  x,  y,  z  to  point 
on  S. 

The  surface  at  Infinity  does  not  c:>ntrlbute  to  the 
value  of  P(xyz),  since  one  of  the  boundary  conditions  Is 
that  the  perturbation  pressure  must  vanish  at  infinity. 
For  the  velocity  potential.  Green's  theorem  gives 


We  shall  consider  both  steady  flow  and  simple  harmonic 
time “dependent  perturbations.  The  steady  flow  case  can 
be  found  by  setting  the  reduced  frequency  k  equal  to 


zero.  The  linearized  relationship  between  the  complex 
amplitude  of  P(xyz)  and  (xyz)  for  simple  harmonic 
motion  is 

(3.5) 

The  inverse  relation  is 

--  ^  ^  (3.6) 

A  sketch  of  the  foil -cavity  surface  S  is  shown  in 
fig.  3-1- 


Fig#  Closed  Surface  of  Cavity  and  Poll 
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As  In  the  linearized  theory  for  fully  wetted  flow,  we 
satisfy  boundary  condition  not  S  but  on  the  projection 
of  S  In  the  x,  y  plane. 

The  boundary  condition  on  the  cavity  surface  Is 
that  the  perturbation  pressure  Is  a  constant.  The 
cavitation  number  0“  is  defined 

T  -  R>r.  -  Pc.  (3.7) 

■ipV- 

pressure  In  the  cavity 
*  free  streemi  pressure 

The  nondlmensional  perturbation  pressure  on  the  cavity 
surface  Is 


Linearized  theory  Is  valid  only  If  cr  1. 

On  the  wetted  surface  of  the  foil  the  boundary  condition  is 

that  there  is  no  flow  through  the  surface. 

The  to'l-cavlty  surface  collapses  to  become  a  region 

In  the  X,  y  plane  as  sketched  In  fig.  3*2.  The  projection 

of  the  foil  surface  Is  denoted  as  S^,  the  projection  of 

the  cavity  surface  as  S  .  When  It  Is  necessary  to 

0 

distinguish  between  the  top  and  bottom  of  S  or  S  ,  we 
will  use  the  notation  S  S  ^  and  S  ~. 
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Using  the  relation  ^(*>5  )  ~  closed 

surface  In  eq.  (3*3)  and  (3*^)  are  rewritten' as  open  sxirface 
Integrals  over  and  S  .  Ihe  pressure  perturbation  at  any 

W  C 

point  xyz  Is  given  by 


5  d 


(3.12) 


idiere 


in 


across  S. 


'  a/  ^ 

across  S  and 


is  the  notation  for  the  Jump 

A  ^ 

A  p  deiiotes  the  Jump  in  p 


(3-13) 

\Sf  /  9f  df 


The  1/K  singularities  are  called  pressure  sources.  It 
can  be  seen  from  the  z  momentum  equation  that  their  strength 
Is  related  to  the  local  curvature  of  the  streamlines.  The 
^ singularities  are  called  pressure  doublets. 

A  pressure  doublet  causes  a  unit  lift  to  act  at  the  point 

^  direction.  For  this  reason,  no 
pressure  doublets  can  appear  off  the  solid  foil  surface. 
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The  expression  for  the  upwash  velocity  perturbation 
caused  by  this  distribution  of  pressure  sources  and 
doublets  Is  found  by  the  application  of  eq.  (3*6)  which 
upon  Interchanging  of  the  order  of  Integration  gives 


c/s  a/' 


7 


7 


The  expression  for  the  upwash  on  the  wetted  foil  surface 
Is 


The  limit  as  of  the  expressions  under  the 

Integral  sign  are  the  Kernel  functions  of  the  problem. 
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Physically  they  represent  the  upnash  at  a  point  on  the 
foil  wetted  surface  caused  by  a  pressure  source  or  pressure 
doublet  of  oscillating  strength  located  at  •  ^e 

analytical  expressions  for  these  kernel  functions  are  evsiluated 
in  Appendix  i.  Starting  from  the  definitions 


.  /  *  iiCx-Cx-g}) 


(3.18) 


The  results  of  appendix  i  give 


(3.19) 


-iUx-3) 

r  /  T  (k{y-3) 

V.  or  k^L  J,  (kii-ni)-L ,  (tlrdiJ  ^ ^ 

//  /  ck(x-3) 

'fi7r 


(3.20) 
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The  kernel  function  Kg  is  used  for  fully  wetted 
flow.  A  discussion  of  the  singularities  of  Kg  Is  given 
In  Chapter  1. 

For  steady  flow  It  becomes 


^  -  /  f,  ^  y-Jt _ 7  (3.21) 

L  -i)j 

The  boundary  condition  on  perturbation  pressure  Is  applied 
by  taking  suitable  llirdta  of  eq.  (3  •12). 


‘  ’’’ fz-  j « ( 


22) 


We  define  a  third  kernel  function 

V  =  y  ^  f .  —  '  -  )  (3.23) 

In  appendix  1  this  Is  shown  to  be 

X-i,  ’  '.2r  (Tf'jc?;  o( -y-s) 
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With  this  'substitution  eq.  (5.22)  simplifies  to 


-  ^(^y) 


(3.24) 


■f  ~^r/f 

The  second  Integral  equation,  for  the  complex  amplitude  of  the 
upwash  distribution  on  the  wetted  foil  surface  becomes 


vfxjy) 


0(3 


(3.25) 


.  // 

ve 

These  two  equations  must  be  solved  for  a  complete  description 


where 


of  the  linearized  supercavltatlng  hydrofoil,  l^ie  unknowns 
are  the  distribution  of  pressure  source  strength 
the  foil  and  cavity  surfaces,  and  the  lift  distribution 
on  the  foil  surface. 

For  steady  flow  at  a  given  cavitation  number  a,  we 
obtain  the  following  equations: 


V(y,y) 


-/ 


■ft?- 


TJ  I 


^  jf  %){cPP>^-y) 

Sw*^t-  ^ 


(3.27) 
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''For  calculation  of  forces  and  isoneit#f  we  assiune 
that  the  unsteady  motion  Is  of  small  amplitude  so  that 
the  cavity  pressure  remains  unchanged  from  Its  steady 
value.  The  following  Integral  equations  for  the  coo^lex 
amplitude  of  the  unsteady  perturbation  then  result: 


(3.28) 


(3.29) 

•a. 


o  = 
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30  Niwerlcal  Method  for  three-dlraenslonal  supercavltating 
Hydrofoil a 

The  coupled  integral  equation  for  the  flow  around  a 
three-dimensional  supercavltating  hydrofoil  were  discussed  In 
section  5.1.  In  the  present  section,  we  shall  describe  a 
numerical  method  similar  to  that  used  for  fully  wetted  foils, 
to  determine  the  lift  distribution  for  steady  and  unsteady 
motion  for  cavities  longer  than  the  chord. 

Equation  (5.26)  and  (5.27)  on  (5.28)  and  (5.29)  will  be 
solved  approximately  by  assuming  a  series  of  functions  for  both 
and  *  The  unknown  coefficients  are  found  by 

satisfying  the  upwash  and  cavity  pressure  boundary  conditions  at 
collocation  points  on  the  foll-cavlty  surface.  For  fully 
wetted  flows  as  discussed  In  Chapter  1,  the  assumed  chordwlse 
functions  for  were  solutions  to  two-dimensional  linearized 

thin  airfoil  problems  Including  the  solution  for  a  flat  plate 
in  steady  flow  with  a  square  root  leading  edge  singularity.  The 
assumed  spauiwlse  functions  were  elliptical  with  going  to 

zero  at  the  tip.  For  an  accurate  determination  of  the  lift 
distribution.  It  Is  Important  to  choose  a  set  of  functions  for 
.hlch  approximate  the  physical  behavior  as  closely  as 
possible.  The  solution  by  linearized  two-dimensional  theory  for 
supercavltating  flow  past  a  flat  plate  has  a  quarter  root  leading 
edge  singularity  in  the  distribution  of  .  We  choose  this 

function  as  one  of  the  modes  ot  ,  in  place  of  the  function 


for  a  fully  wetted  flat  plate.  'Hie  integrations  of  these 
assumed  modes  time  the  kernel  function  relating  upwash  to 
are  preformed  numerically. 


Because  of  the  very  simple  fona  of  the  kernel  function 
relating  pressure  in  the  cavity  to  the  pressure  source  strength 
distribution,  we  chose  functions  for  which  the 

Integrals  In  eq.  (3.27)  and  (3.29)  can  be  evaluated  analytically. 

Since  the  distribution  Is  affected  by  the  ^ 

distribution  only  through  the  coupled  integral  equations,  we 
should  be  able  to  obtain  an  accurate  prediction  of  lift  and 
moment  with  only  an  approximate  pressure  source  strength 
distribution.  Rather  than  select  as  streamwise  modes  of 
the  distributions  of  from  linearized  two-dimensional  solutions 

for  supercavitatlng  flows,  we  use  a  simpler  set  of  functions 
which  replaces  the  actual  singular:! ties  of  these  solutions 
with  simple  delta  functions  and  gives  a  piecewise  linear 


approximation  to 


(g). 


The  two-dimensional  linearized  solution  for  a  supercavitatlng 
flat  plate  has  a  5/4  root  leading  edge  singularity.  The  I 

correct  definition  of  Integrals  over  this  singularity  is  the  I 

Kadamar  "finite  part"  (ref,  (11)).  Since  the  velocities  amd 
pressures  produced  by  this  non-lntegrable  singularity  have  only 
a  quarter  root  leading  edge  singularity,  we  feel  that  the 
boundary  conditions  can  be  satisfied  and  an  accurate  prediction 
of  obtained  by  approximating  the  effect  of  this  singularity 


<lp 


with  a  piecewise  linear  l\lnctlon  plus  a  delta  function. 
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Both  the  supercavitating  flat  plate  and  symmetric  wedge 
solutions  have  a  square  root  singularity  in  at  the 

trailing  edge.  We  hope  to  approximate  the  effect  of  these 
singularities  of  oh  the  lift  distribution  by  locating 

delta  function  of  arbitrary  strength  at  the  leading  and 
trailing  edges.  Although  this  is  admittedly  a  very  crude 
procedure,  it  is  felt  that  it  would  be  sufficiently  accurate 
at  the  present  exploratory  stage. 

In  the  two-dimensional  non-linear  theories  of  finite 
cavity  flows,  various  au?tif ices  ha.e  been  used  to  account  for 
the  trailing  edge  of  the  cavity.  The  constant  pressure 
condition  on  the  cavity  prevents  the  occurence  of  a  stagnation 
point  with  smooth  flow  off  the  trailing  edge.  • 

Two  models  which  have  been  used  for  finite  cavities  are 
the  re-entrant  jet  model  and  the  modified  Riabouchinsky  model 
which  terminates  the  cavity  with  a  fictitious  vertical  flat  plate. 


Modified  Riabouchinsky  model 
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Linearized  theory  derived  from  both  of  these  models 
replaces  the  trailing  edge  of  the  cavity  with  a  singularity 
and  imposes  a  cavity  closure  condition. 

Another  non-linear  model  which  has  been  studied  is 
the  smooth  wake  termination  model. 


The  linearized  theory  based  on  this  model  does  not  require 
cavity  closure,  but  only  that  the  vertical  velocity  is 
continuous  in  the  wake. 

Pabula  (ref.  (55))  has  compared  the  lift  and  moment 
computed  by  these  linearized  theories  derived  from  these  two 
types  of  cavity  models.  For  cavities  at  least  twice  as  long 
as  the  chord,  that  is  for  9c<  1>  the  numerical  results 

differ  by  less  that  2jt. 

In  terras  of  pressure  source  strength,  the  cavity  closure 
condition  would  require 

/ 

The  smooth  wake  termination  would  require 


=  <D 
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Either  of  these  conditions  could  be  incorporated  into 
the  numerical  method  by  adding  a  cavity  termination  constraint. 
Further  research  should  be  done  to  determine  the  effect  of 
closure  conditions  at  the  trailing  edge.  It  is  not  expected 
that  the  actual  conditions  used  will  effect  the  prediction 
of  lift  and  moment  for  cavities  longer  than  the  chord. 

In  the  present  method,  we  terminate  the  cavity  by 
requiring  the  pressure  source  strength  to  go  to  zero  beyond 
the  cavity  trailing  edge. 

Further  refinements  could  be  made  to  account  for  the 
actual  singularities  by  choosing  as  the  chordwise  modes  of 
the  solutions  from  linearized  two -‘dimensional  theory. 

The  integrations  would  then  be  performed  by  by  numerical 
quadrature.  The  present  discussion  merely  suggests  a  simple 
numerical  technique  for  three-dimensional  supercavitating  foils. 
The  success  of  the  approximations  can  be  Judged  by  the  accuracy 
of  the  prediction  of  lift  and  moment  on  the  foil. 
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3*  3a  Numerical  Calculations  for  a  Steady  Synvnetrlc  Wedge 

We  shall  begin  our  discussion  of  numerical  techniques  for 
three-dimensional  supercavi bating  flows  by  considering  the  epe^al 
case  Of  a  symmetric  wedge. 

For  this  special  case,  the  integral  equations  (3*26)  and 
(3.27),  uncouple.  Since  Is  lero  for  synmetric  wedge 

flew,  equation  (3*26)  becomes  simply 


C  9% 


The  solution  for  a  wedge  of  half  angle  <k,  is 


(5.30) 


(3.31) 


The  second  integral  equation,  (3*27)  becomes 


(3.32) 


The  only  unknown  is  the  distribution  of  ®c' 

One  additional  problem  in  cavity  flows  is  that  the  shape  and 
length  of  the  cavity  is  not  known  In  advance.  To  proceed  with  a 
numerical  method,  however,  we  must  make  some  assumptions  about  cavity 
length  and  shape.  We  then  test  our  numerical  results  to  determine 
Che  sensitivity  to  these  assumptions.  For  a  valid  linearized 
theory  we  require  the  cavir.ation  number  to  be  smsdl.  In  the  exact 


oi 


limiarized  solution  for  two  (Unensional  cavity  flows,  long  cavities 
are  consistent  with  small  cavitation  nunbers.  Thus,  any  effects 
due  to  cavity  length  and  end  conditions  should  produce  minor 
effects  on  the  foil  itself. 

We  assume  that  the  cavity  surface  is  a  rectangle  equal  in 
span  to  the  foil  and  of  some  length  X^.  For  a  given  foil  and  a 
given  cavitation  number  w  would  test  the  numerical  results 

for  sensitivity  to  the  assumed  value  of  X^.  Fig.  shows  the 
configuration  for  the  linearized  solution  of  a  siq>ercavitating 
symmetric  wedge  of  constant  chord  equal  to  two,  and  of  aspect  ratio 
S.  In  the  linearized  problem,  both  the  Voil  and  cavity  surfaces 
are  located  in  the  x,y  plane  (z  »  0). 


Fig.  Model  for  a  supercavitating  wedge 

For  this  problem,  the  source  distribution  ®  symmetric 

function  about  the  axis. 

As  in  the  numerical  solution  for  fully  wetted  three  dimensional 
flows,  we  pick  a  set  of  n  functions  for  the  source  distribution  with 
unknown  coefficients.  These  functions  will  be  symmetric  in  . 

The  n  coefficients  are  found  by  satisfying  eq.  (507)  at  collocation 
points  on  the  cavity  surface. 
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Ve  define  a  set  of  spamfiee  functione 
the  following  properties: 


j-J.- 

T'  ^ 

7 

^  V,  ^  -Vi. 

%  ^ 

-^Z 

^  7*  ^  t; 

^  K  Z 

w  - 1 

I'l*  " 

% 

tz 

7,  *tx 

/ 

l-t*  ^ 

7z  “^7 

(3.38) 

The  fxmction  sketched  in  Fig*  5.4 


-"1,  "1,  "t.  V  \ 

Pig.  5*4  rtie  function  '*1 

The  rectangular  cavity  areSf  &.»  is  divided  into  km  rectangles 
of  width  length  •  Ve  let  I  be  the  niaber  of 

rectangles  in  the  3  direction  and  SET  be  the  niaber  of  rectangles 
in  the  direction.  Then 

A-j  -  i  Xe  -  0  /  X 

^  7'  (3.39) 

The  following  set  of  functione  for  are  chosen: 

'3f  ' 

J  ^  r 


(3.'*0) 


c 


where 


Ihe  f line t Ions 


ict)  -  Tct-O  -»  ^2 


(J.'tl) 


are  sketched  In  Pig,  3.5, 


Pig.  3.5 


The  Set  of  Functions  -^-  c-i) 


S3 
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The  function  has  the  following  properties: 

1.  Unit  delta  function  behavior  at  the  trailing  edge, 
where  the  two  dimensional  linearized  solution  has  a 
square  root  singularity. 

2.  Piecewise  linear  and  continuous  variation  in  ^  from 
the  trailing  edge  to  the  end  of  the  cavity. 

For  a  typicail  set  of  O.:  >3  ,  the  function  is  sketched  in 


The  assumed  set  of  functions  for  \0xjrv;  has  the 
following  properties, 

1.  Symmetric  in  ,  constant  in  each  rectangle. 

2,  Symmetric  delta  function  rows  of  unknown  strength  at 
the  trailing  edge  in  each  rectangle. 

Piecewise  linear  and  continuous  in  3  . 


The  total  number  of  unknown  coefficients  is  (J  T 
which  is  also  the  number  of  collocation  points  required.  We  put 
two  control  points  in  each  rectangle  closest  to  the  trailing  edge 
to  allow  for  accurate  determination  of  the  delta  function  strength 
and  to  satisfy  boundary  conditions  close  to  the  trailing  edge.  The 
remaining  control  points  are  located  at  the  '  coordinates  in  the 
center  of  the  reotanglesu 


(See  Pig.  J.7  ) 


5,7  Control  point  locations  on  for  a  syimetrlc 
wedge  solution 

These  siisple  assumed  distributions  allow  the  integral 
in  (5.37)  to  be  evaluated  analytically*  For  the  n  collocation 
pointSi  we  thus  obtain: 


P 


Kr 


(3.42) 


where 


Wt  -  ± 

air 


coordinates  of  the  kth  control  point  and  Pj^ 
is  the  pressux^  at  the  Kth  control  point  caused  by  the  n^eji 
assumed  ^  mode. 

Solution  of  equation  (3*42)  gives  the  s  »  the  coefficients  of 

the  source  distribution  modes.  The  pressure  at  any  point  on  the 
wedge  may  be  found  from 

^  IrCx,^)  <^n 

where  Pn  (xy)  is  the  pressure  at  a  point  x,y  due  to  the  nth  source 
mode,  and  W  (xy)  is  the  pressure  at  a  point  x,y  due  to  the  leading 
edge  delta  function  source  of  strength  2o{, 
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J.Jb  Numerical  Reeulte  for  Sytaoietrlc  Wedge 

The  numerical  method  has  been  programmed  for  the  IBM  7090* 
Ve  will  present  a  few  numerical  results,  k  useful  comparison  can 
be  made  between  numerical  results  obtained  for  pressure  and 
source  strength  distributions  on  the  center  section  of  a  large 
aspect  ratio  wedge  and  results  of  two  dimensional  linearized 
theory  for  cavltatlng  wedges. 

Two  dimensional  linearized  theory  for  supercavltatlng 
wedges  at  zero  cavitation  number  gives  the  following  results 
(see  Guerst  ref.  Q2) ) 

Hie  pressure  distribution  along  a  symmetric  wedge  of 
semichord  unity  Is: 


(5.^5) 


where  x  >  >1  Is  the  leading  edge 
X  >  -t-l  Is  the  trailing  edge 
and  Is  the  wedge  half  angle  In  radians. 

The  pressure  source  distribution  Is  zero  on  the  wedge 
except  for  a  delta  function  of  strength  2^  at  the  leading  edge. 
In  the  cavity  the  distribution  of  soiurce  strength  Is 


-- 

9f  ' 


(5.44) 


Results  were  computed  for  an  aspect  ratio  6  wedge  with  half 
angle  o<  •  .1  radians.  Hie  calculations  were  done  for  cavity 
lengths,  of  11,  7,  ^  and  3  at  a  cavitation  nunber  of  zero. 


0  0,1 


C 

FIG.  3.8  COMPARISON  OF  CHORD^^ISE  PRESSURE 

DISTRIBUTION  ON  A  SYMMETRIC  WEDGE  BY 
TWO-DIMENSIONAL  THEORY  AND  THREE-DIMENSIONAL 
NUMERICAL  RESULTS  ON  THE  CENTER  SECTION 
OF  AN  ASPECT  RATIO  6  WEDGE  a‘0  a«O.I 


3  5  7  9  II 


FIG. 3.9  COMPARISON  OF  PRESSURE  SOURCE  DISTRIBUTION 
IN  THE  CAVITY  BEHIND  A  SYMMETRIC  WEDGE  BY 
TWO-DIMENSIONAL  THEORY  AND  THREE-DIMENSIONAL 
NUMERICAL  RESULTS  a  =  0  o  =  0.1 
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Fig.  8  shows  a  comparison  between  the  pressure  distribution  on 
the  center  section  and  the  two  dlnensionail  theoretical  result. 

The  numerically  calculated  cuz*ve8  for  various  cavity  lengths  fall 
on  top  of  one  another,  and  thus  are  very  insensitive  to  assumed 
cavity  length. 

Pig.  3*  9  gives  a  comparison  between  the  source  strength 
on  the  cavity  surface  of  the  center  section  as  calculated 
numerically  for  various  assumed  cavity  lengths  and  the  two 
dimensional  result. 

3.3c  Numerical  Method  for  a  Supercavitatlng  Flat  Plate 

The  numerical  approach  to  the  supercavitatlng  flat  plate 
is  similar  to  that  described  for  the  symmetric  wedge.  In  this 
case,  however,  the  integral  equations  do  not  uncouple  and  must  be 
solved  simultaneously. 

The  cavity  shape  is  again  taken  as  a  rectangle  on  the  Z  •  0 
plane  of  some  length  and  of  span  equal  to  that  of  the  foil. 

In  this  problem,  there  are  two  unknown  distribution  of 
singularities,  a  pressure  doublet  distribution  on  the  foil  only 
and  a  pressure  source  distribution  on  the  foil  and  cavity  surfaces. 
For  each  of  these  distributions,  we  pick  a  set  of  functions  with 
unknown  coefficients.  The  boundary  conditions  on  both  upwash  and 
cavity  pressure  are  then  satisfied  at  an  equivalent  number  of 
appropriate  control  points  on  foil  and  cavity  surface.  Results  will 
yield  lift  distribution  on  the  foil  and  pressure  source  distribution 
of  foil  and  cavity. 
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The  assumed  series  for  was  that  used  for  fully 

wetted  flow.  Itils  series  allows  the  possibility  of  a  square  root 
singularity  at  the  leading  edge.  Two  dimensional  solutions  for 
supercavitating  flat  plates  show  a  quarter  root  singularity  at  the 
leading  edge.  From  the  numerical  standpoint,  it  is  a  simple 
matter  to  substitute  this  function  for  the  square  root  singularity. 
For  steady  flow  solutions,  we  consider  only  functions  which  are 
symmetric  in  y  oscillation  about  a  steady  cavity  solution, 

we  consider  both  symmetric  and  antisymmetric  solutions.  Itie  form 
of  ,  as  in  chapter  1,  is  taken  to  be 


V  tA 


where  -f  ~  O 


'"iLI 


/'ll' ) 


'y. 


(J.'ts) 


(3.46) 


the  symmetric  functions  are 


(3.47) 


'i I 

the  antisymmetric  functions  are 


(3.48) 


z. 


l  r  ■*  I 


(3.49) 
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The  Integrals  in  equation  (3.?3)  and  (3*35)  involving 
£i.p(3  are  treated  numerically  by  the  same  technique 

used  in  the  numerical  method  of  chapter  1  for  fully  netted  flow. 

The  assumed  set  of  functions  for  the  pressure  source 
distribution  on  the  foil  and  cavity  surfaces  is 

aimilar  to  that  chosen  for  the  symmetric  wedge  solution.  Both  the 
foil  and  the  cavity  surfaces  are  divided  into  rectangular  areas 
of  width  and  length  on  the  foil  and  on  the  cavity. 

and  ax*e  not  necessarily  equal  in  length.  The 

form  of  <¥M  is  chosen  as  follows: 

For  steady  flow  cavity  solutions  ,dr  for'^a  symmetric 
oscillation  about  a  steady  flow  solution  is  symmetric 

in  /y  . 


(3.50) 


where  h  (  ^  a  p  )  defined  in  equation  (3.38). 

For  antisymmetric  oscillations  about  a  steady  cavity  flow, 
(^  ^)\  antisymmetric  in  .  We  therefore  set 


=■  ^  (3.51) 


where 
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Since  the  two  dineneional  linearized  solution  for  a  super- 
cavitating  flat  plate  has  singularities  in  {^ ) 
both  the  leading  and  trailing  edges,  the  function  f (3  )  1b 
chosen  to  have  the  following  properties: 

1.  unit  delta  function  at  the  leading  edge 

2.  arbitrary  strength  delta  function  at  trailing  edge 
3*  piecewise  linear  and  continuous  on  both  foil  and 

cavity  surfaces,  hay  be  discontinuous  at  the 
trailing  edge. 

If  the  number  of  rectangles  in  the  3  direction  is  l^ 
on  the  foil  and  on  the  cavity,  the.  total  nmber  of 

inodes  in  f(s)  (if  2)  v  (l^  +  2),  These  modes  of 
f(3)  sketched  in  Fig, 


.ju)  -  Sc3i-|)  +• 

I_^lc  +  4 
4- 


di. 


)  (3.52) 


For  a  set  of  d,'B#^a  typical  fC*^)  iB  sketched  in  Fig.  3*10 . 

.A  — ' 


t 


Fig,  5.10 


Typical  f(3)  for  supercavitating  plate. 


The  Integrals  of 
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^  in  equations  (3*33)  to  (3.35) 
can  be  done  analytically  for  this  set  of  modep. 

The  total  number  of  rectangular  areas  on  S  amd  S  is 

Vrf  W 

2N  *  (l^  -f  I  ).  The  total  number  of  unknown  coefficients  for  the 

X  c 

presijure  source  modes  Is  N  ♦  (I^  +  2  +  +  2).  For  nianerlcal 

simplicity,  we  choose  N  spanwise  modes  end  2  chordwise  modes 

for  on  the  foil  surface.  We  then  satisfy  both  the  upwash 

condition  and  the  cavity  pressure  condition  at  the  same  control  points 
on  the  foil.  We  need  N  ♦  (I^  2)  control  points  on  the  foil  for 

either  the  symmetric  or  anti synsse trie  case.  The  condition  on  cavity 
pressure  behind  the  foil  on  is  satisfied  at  N  ♦  (I^  +2)  control 
points.  This  collocation  gives  a  complete  set  of  linear  algebraic 
equations  whose  unknowns  are  the  coefficients  of  the  modes  and 

the 

The  set  of  control  points  u^ed  In  the  numerical  method  is 
sketched  In  Pig.  3.  12  .  The  philosophy  on  control  point  location 

is  to  select  an  even  distribution  with  special  “raphasls  on  sensitive 
areas  such  as  leading  or  trailing  edges.  For  either  symmetric  or 
antisymmetric  cases,  control  points  are  located  on  bnly  half  of 


Fig.  3.12 


Configuration  for  numerical  solution  for  super- 
cavltatlng  flat  plate  showing  control  point 
locations, 
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In  matrix  notation,  equations  (5* 35)  and  (5.54)  on 
(5.55)  and  (5.56)  take  the  following  form  after  the  Integrals 
have  been  evaluated  either  numerically  or  analytically  for  the 
assumed  set  of  and  modes. 


DW 

DWW 

0 

-cr 

2. 

— 

iC»-. 

P 

~<r 

a 

0 

> 

k 

(^.55) 


where  Is  the  coefficient  of  the  n^j  mode  on 

the  foil  and  a^^  Is  the  coefficient  of  the  n^i^  mode. 

The  matrix  has  been  partitioned  to  show  the  relation¬ 
ship  to  equations  (505)  and  (5.54)  or  (5.55)  dJiii  (5*56).  For 
the  latter  set  of  equations  (l.e..  In  unsteady  flow)  this  matrix 
set  Is  complex.  The  elements  of  the  psirtltlon  DW  are  the 
upwash  at  the  kth  control  point  from  the  nWj  mode. 

For  a  fully  wetted  flow,  an  Inversion  of  this  matrix  alone  gives 
the  complete  solution  for  the  lift  distribution. 


■'  1 

Vk 

DW 

b.  ^ 

«  — 

(3.54) 
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The  elements  of  the  IU)V  matrix  are  found  from  the  integral 
on  in  eq[uat5on  (3.33)  or  (3*35) 


DDw  =  "  -V  I  ^ 


(3.55) 


Only  those  3  points  directly  ahead  of  the  control 
point  contribute  to  the  elements  of  WU»  Ttie  elements  of  the 
DPM  are  just  the  modes  evaluated  at  the  control  points  of  the 


foil. 


DPM 


(3.56) 


where  denotes  the  nt^  pressure  mode. 

Hie  elements  of  the  P  matrix  are  the  pressures  at  the 
control  points  on  both  the  cavity  and  foil  surface  due  to  the 
various  assumed  modes  ^ 


of 


F  =  - 


^ir' 


cJ'^  (/ 


(5.57) 


where  source  node  and  >5^^  is  the 

appropriate  area  for  the  n^  source  mode. 

The  two  matrices  with  zero  elements  appear  in  the  complete 
matrix  for  two  reasons:  first  the  pressure  source  distribution  on 
the  cavity  surface  does  not  affect  the  upwash  .  on  the  foil 
surface,  second  the  function  is  zero  on  the  cavity 


surface. 


The  matrix  in  eq.  (3.33)  is  inverted  to  give  the 
coefficients  b^'s  and  a^»B  for  a  given  cavitation  number 


96 


and  upwash  distribution  v(xy)  on  the  wetted  surface  of  the 
foil.  Hence { 


(3.58) 


where  the  fc]  matrix  is  the  inverse  of  the  square  matrix  in 

(3.53}.  For  a  flat  plate^  the  upwash,  v(x,y)  is  equal  to 

-  ^  •  Ve  rewrite  equation  (3*58)  to  show  th«  dependence  of 

the  a  *  s  and  b  •  s  on  o(  and  0“ 
n  n 


(3.59) 


(3.60) 


Cl  jy  denote  the  partitions  of  the  C  matrix.  For  a 

given  assumed  cavity  length,  the  coefficients  are  linear 
in  ^  and  ^  with  coefficients  proportional  to  the  sum  of 
row  elements  of  [C^  • 


htr.  ^  y-  •  (T 

'<y<  'C~ 


(3.61) 
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r  f  f  - 
i.„x  '  f 

ior  - 

Ttie  lift  and  aonent  on  the  foil  aay  be  obtained  froa 
a  linear  operation  on  the  b^*8  aa  in  the  case  for  a  fully 
wetted  surface. 


r  L  J  C 


(5.62) 


J 


u 


j 


(3.65) 


where  and  C_  are.  row  na trices;  C,  #  C_  are  the 

“ti  “n 

contributions  to  the  lift  and  moment  coefficients  of  the  n^ 

mode.  The  lift  and  moment  coefficients  are  then  linear 
in  both  o<  and  cr*  • 

Q  -  ■  ‘T 


where 


(3. 64 


C4/ 

• 

1 

j’ 

L 

r>» 

X 

S 

-i. 

buflj 


(3.65 
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For  any  given  foil,  the  coefflcente  C,  ,  and  c, 

1j  cA  Li  r 

may  be  calculated  for  various  values  of  the  assumed  cavity 
length  Xg  to  determine  the  sensitivity  of  the  predicted  loads 
to  cavity  length.  As  in  the  case  of  the  non-lifting  wedge, 
the  results  for  the  lifting  flows  considered  below  were, 
fortunately,  found  to  be  very  insensitive  to  cavity  length, 
provided  is  longer  than  about  twice  the  foil  chord. 

Results  of  two  dimensional  linearized  theory  and  the  experiments 
of  Kermeen  (ref  (20))  indicate  that  the  cavity  is  long  compared  to 
the  chord  for  ^  j_  .  This  must  be  kept  in  mind  when 

applying  the  numerical  method. 

The  numerical  method  described  in  this  section  has 
been  programmed  for  the  IBM  7090  for  both  steady  and  unsteady 
flow.  Ve  will  present  a  few  numerical  results  and  comparisons 
with  both  theory  and  available  experiments. 


5.  >d  Comparison  with  Steady  Two  Dimensional  Linearized  Theory 

A  useful  comparison  can  be  made  between  two  dimensional 
linearized  theory  and  some  numerical  results  for  the  center 
section  of  a  large  aspect  ratio  supercavitating  foil,  in  this 
case  an  aspect  ratio  6  foil. 

Two  dimensional  linearized  theory  gives  the  following 
results  for  lift  distribution  on  a  steady  supercavitating  flat 
plate  at  a  cavitation  number  of  zero  (Ref.  ^2). 

/  /  ( 

Ap('v)  =  o<  f  -  ^  ^  / 


(3.66) 
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where  ^  Is  the  angle  of  attack 

X  >  -X  is  the  leading  edge 
and  X  >  4-1  Is  the  trailing  edge 

The  expression  for  the  pressure  source  strength  on 
this  foil  Isi 


(5.67) 


where  the  "finite  part"  must  be  taken  of  Integrals  over 
the  5/i|  root  singularity. 

The  source  strength  on  the  cavity  Is 


These  expressions  are  plotted  In  Pig.  5.15  and  5. lU  for 

^  -  .1. 

Calculations  using  the  present  method  were  made  for  a 
rectangular  flat-plate  foil  of  aspect  ratio  6  at  ^  -  0. 

Results  for  the  chordwlse  lift  and  pressure  source  distribution 
at  the  center  section  are  Included  In  Figures  3.13  and  3.14 
For  these  calculations,  a  cavity  length  equal  to  3  times 
the  chord  was  chosen. 

Calculations  were  done  with  both  a  square  root  and 
quarter  root  leading  edge  singularity  In  the  assumed  form 
of  •  Fig.  3.13  fthd  3.14  also  compare  the  results 

from  these  two  methods*  The  difference-in  predicted  lift 


FIG.  3.13  NUMERICAL  RESULTS  FOR  THE  PRESSURE  SOURCE 

DISTRIBUTION  ON  THE  CENTER  SECTION  OF  AN  ASPECT 
RATIO  6  SUPERCAVITATING  FLAT  PLATE ,  COMPARISON 
WITH  TWO-DIMENSIONAL  LINEAR  THEORY, 
a  =  0.1  ,  (T  =  0,  X£  =  5 
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FIG  3  14  NUMERICAL  RESULTS  FOR  THE  LIFT  DISTRIBUTION  ON  THE  CENTER 
SECTION  OF  AN  ASPECT  RATIO  6  FLAT  PLATE  ,  COMPARISON 
WITH  TWO -DIMENSIONAL  THEORY  a=O.I,  o-=0,Xg  =  5.0 


TABLE  3.1  NUMERICAL  RESULTS  FOR  THE  CENTER  SECTION  OF 
AR  =  6  FLAT  PLATE  ,  COMPARISON*  WITH  2- D  THEORY 


2-0 

NUMERICAL 
1/4  root  sing. 

NUMERICAL 
1/2  root  sing. 

c. 

•-a 

1  .5707 

1 . 52 

1 .6055 

^CP 

0.32 

0.  29 

0.275 
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coefriclent  is  .al)out  The  variation  in  the  center  of 

pressure  location  iS  also  about  5^.  The'  prediction  of 


source  strength  ie  very  orude^  but  seeMS  adequate  for  load  calculation! 
Two  dimensional  linearized  (Ref.  >2)  theory  gives 


the  following  results  for  lift  and  moment  coefficient  and  cavity 


length. 


y.  =  ) 


(3.69) 


^ L  ~  '^iJT  ^ ^  )  (  I  ^ 

ie  a 


(3.70) 


C>y,  --  C”* 


^  CoS  C 0 


8  (  CcsYc  ( if^CosT,)"^. 


(>.71) 


^  ^ 


(5.72) 


where 


t  1 


a  a  parameter 


The  expressions  for  C_  and  C  for  o<  »  ,l  are  plotted  vs,  'T 

L  m 


In  Fig.  5,15  The  nvanerical  results  obtained  for  the 


and  C  on  the  center  section  of  the  aspect  ratio  6  foil  are 


also  plotted  vs.  cr  .  As  discussed  in  section  3,3c.  the 


numerical  technique  gives  a  linear  relation  for  and 


(3.73) 


Cm  =  .■<<  ^  Cm 


(3.7“) 
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The  non-linear  character  of  the  C*  and  C  va  curves 
predicted  by  two  dimensional  theory  cosies  from  the 
dependence  of  cavity  length  on  <7*  • 

The  three  dimensional  results  have  shown  a  very  small 
variation  of  these  coefficients  wiUi  assumed  cavity  length  • 
Table  5.2  gives  Ct  and  C,  _  tabulated  vs. 

assumed  cavity  length 


FIG.3I5  LIFT  AND  MOMENT  COEFFICENTS  VS  cr,  (a  =  0  1) 

NUMERICAL  RESULTS  FOR  THE  CENTER  SECTION  OF 
AN  AR  =  6  FOIL;  TWO-DIMENSIONAL  THEORY 


TABLE  3.2  NUMERICAL  RESULTS ,  TABULATED  FOR  ASSUMED 
CAVITY  length 


o 

r~ 

n 

mam 

mm 

5.0 

1 .5131 

.4019 

-  44159 

■a 

7.0 

1 . 5097 

.3783 

-  44058 

mm 

9  0 

1  5000 

.3629 

*  4300 

-.1337 

u,  i- _ i. 


0  +1 


■» 
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3.3e  Coroparison  with  Other  Three -Dimenalonal  Theories  and 


5 


Experiment 


Johnson  (ref.  (21),  (22))  has  developed  a  method  for 
calculating  lift  and  drag  on  a  supercavltatlng  hydrofoil  at  a 
cavitation  number  of  zero.  This  theory  Is  based  on  the 
assumption  tha:  the  "Influence  of  finite  span  on  the  two- 
dimensional  lift  coefficient  Is  due  to  the  effects  of  the 
trailing  vortlclty".  For  fully  wetted  flow,  this  assumption 
forms  the  basis  of  lifting  line  theory  and  Its  various  refine¬ 
ments. 

The  results  of  this  assiamptlon  In  fully  wetted  flow  give 
the  following  form  for  the  lift  coefficient,  (see;  e.g.,  Jones 
Ref.  (24)). 


Cl.  *  Zfr  ( i') 


(3.74) 

where  Is  the  Induced  angle  due  to  the  trailing  vortlclty. 
Johnson  assumes  that  for  a  supercavltatlng  hydrofoil,  the  lift 
coefficient  takes  the  same  form  with  2Tr  ,  the  lift  curve  slope 
for  two-dimensional  fully  wetted  flow  replaced  by  Tr/2* 
lift  curve  slope  for  supercavltatlng  two-dimensional  flow  at 
a  cavitation  number  of  zero. 


(J.75) 


evaluated  from  lifting  line  theory  for  elliptical  loading  Is 


Q/ .  _ 


(3.76) 


Glauerts  modification  (Ref.  (23) )  to  is  used  for 

rectangular  planfoms. 
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(//r) 

if^  (?-77) 

where  is  a  planform  correction.  A  correction  for  lift  due 
to  crossflow  around  the  tips  is  also  included 


/ 

/ 


-  '^5  'K 


(3.78) 


The  second  theory  for  steady  supercavitating  foils  was 
developed  by  Cumberbatch  Ref.  (23).  'Hiis  theory  is  valid  for 
large  aspect  ratio  foils  at  cavitation  nunbers  for  which  the 
cavity  length  is  of  the  order  of  the  chord,  (-I*  Xg‘  1.3). 

The  model  developed  for  this  theory  is  a  large  aspect  ratio  foil 
with  two-dimensional  supercavitating  flo».  over  the  foil  except 
at  the  tips.  At  the  tips,  the  wing  tip  vortices  cause 
conical  wing  tip  cavities.  These  wing  tip  cavities  are  treated 
by  slender-body  theory.  The  effect  of  the  wing  tip  vortices 
or.  the  flow  away  from  the  tips  is  handled  by  an  induced  angle 
correction. 


Ci_  ^  C.  cx'-Vf) 

2  C 


(3.79) 


where  is  the  induced  angle  created  by  the  two  trailing 

cavltating  vortices.  W  the  two-dimensional  lift 

coefficient  at  cavitation  number  cr  and  local  angle  of  attack 


e 
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A  few  steady  flow  calculations  were  perforated  to 
compare  with  the  experiments  of  Scheibe  and  Wetsel  (Ref.  19) 
and  with  those  of  Kermeen  (Ref.  20).  The  first  set  of 
experiments  measured  lift  coefficients  on  flat  plate 
supercavitating  foils  of  aspect  ratio  2.3#  ^*0,  and  6.0. 

The  cavities  were  artifically  created  by  injecting  air  on 
to  the  upper  surface  of  the  foil.  'Ihe  range  of  CT  obtained 
in  these  tests  was  from  about  .03  to  .3.  In  the  second 
set  of  experiments^  lift  coefficients  were  measured  for 
supercavitating  flat  plate  foils  of  aspect  ratio  4,2,1 
and  l/2.  The  cavities  on  the  foils  were  obtained  naturally 
by  pumping  down  the  ambient  tunnel  pressure.  Two  values 
of  0"  were  obtained  for  each  test,  CTy^  calculated  from  the 
veipor  pz*essure  of  water  at  the  test  conditions  and  the 
actual  measured  in  the  test. 

Photographs  included  in  both  of  these  references 
reinforce  our  choice  of  a  model  for  the  cavity.  The  cavities 
are  rectangular  in  shape,  bounded  by  the  cavities  created 
by  the  trailing  vortices.  The  trailing  edge  of  the  cavity 
is  not  sharply  defined,  but  disappears  in  a  region  of 
turbulence  and  bubbles.  Our  neglect  of  a  strong  trailing 
edge  singularity  seems  Justified  by  these  photographs. 

Experimental  results  of  cavity  length  v.s.  (^) 
are  shovm  in  Fig.  3.16  •  If  the  cavity  length  is  about 

three  times  that  of  the  chord,  our  nuaerical  method  is 


valid 


lo6 


Pig.  3A7  Bhows.  the  results  of  nuserlcal  calculations 
for  an  aspect  ratio  6  flat  plate  supercavltatlng  foll»  as 
compared  to  experiment.  Two  cavity  lengths  were  assumed^ 
as  cavity  of  three  times  the  chord,  and  a  cavity  of  five 
times  the  chord.  'Ate  results  of  the  calculation  for 

m  .148  show  the  insensitivity  of  the  load  prediction 
to  assumed  cavity  length.  The  lift  coefficients  for 
cr  •  0  as  predicted  by  the  method  outlined  by  Johnson  are 
also  shown  on  this  figure.  Since  this  method  is  based 
on  linearized  theory,  both  and  ^  should  be  small; 
the  disagreement  at  high  and  ^  is  not  unexpected. 

Calculations  were  done  for  an  aspect  ratio  4 
supercavltatlng  foil  at  an  angle  of  attack  of  10^ 

(  «»  .175  rad.)  with  assumed  cavity  length  three  times 

the  chord.  Pig.  3.18  shows  lift  coefficient  vs 
cavitation  number  for  the  numerical  results  obtained  by 
this  method  as  compared  to  the  theory  of  Cinberbatch 
(ref. (23))  and  with  the  experiments  of  Ref,  (19)  and  (20). 
The  T  m  0  point  as  predicted  by  the  method  of  Johnson  is 
also  shovm.  The  numerical  method  is  in  good  agreement  with 
the  available  theories,  the  agreement  with  experiment  is 
13  rather  poor. 


CAVITY  LENGTH/CHORO  LENGTH 


CAVITATION  NUMBER  ^ 

ANGLE  OF  ATTACK  IN  RADIANS 


FIG.  3,16  EXPERIMENTAL  RESULTS  FOR  CAVITY  LENGTH 
VS  cr/a 
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FIG  3  17  COMPARISON  BETWEEN  NUMERICAL  RESULTS  AND  EXPERIMENT 


FOR  AN  ASPECT  RATO  6  FOIL  LIFT  COEFFICENT  VS  CAVITATION 
NUMBER  (T 


0  0  1  0.2 


CAVITATION  NUMBER  -  cr 

FIG  3  18  LIFT  COEFFICENT  VS  CAVITATION  NUMBER  FOR  AN  ASPECT 
RATIO  4  FOIL  (a=  I0«),  COMBARlSiON  WITH  EXPERIMENT  AND 
OTHER  THEORIES 
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50f  Lift  and  Moment  on  an  Oscillating  Supercavitating 
Hydrofoil  of  Finite  Span 

The  lift  and  moaent  coefficents  for  an  oscillating 
supercavitating  hydrofoil  were  calculated  using  the 
numerical  method  described  previously.  For  these  cases,  the 
boundary  conditions  on  the  complex  amplitude  of  upwash  on 
the  wetted  surface  was  satisfied  at  13  control  points,  the 
boundai*y  condition  on  Uie  pressure  In  the  cavity  was  satisfied 
at  33  control  points.  Results  for  assumed  cavity  lengths  of 
two  and  three  times  the  chord  differed  by  about  2%.  Compu> 
tatlon  time  for  a  case  with  two  assumed  cavity  lengths  was 
under  three  minutes  on  the  IBM7090. 

The  unsteady  aerodynamic  coefflcents  for  the  center 
section  of  an  aspect  ratio  six  foil  In  heave  and  pitch 
about  the  leading  edge  are  compared  with  the  two-dimensional 
coefflcents  as  predicted  by  Woods  (ref.  29).  These  coefflcents 
are  plotted  vs  reduced  frequency  k  in  Fig. 3.19  to  3.22 
by  use  of  the  following  definitions; 

lift  due  to  heave 

lift  due  to  pitch  about  L.E. 

Z'-  *  ha 

moment  about  L.E.  due  to  heave 
moment  about  L.E.  due  to  pitch 


cxl 


'’hr  '^i 


Cjl  i"  1  Cy  • 

"txr  "oM 
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The  total  value  of  the  unsteady  aerodynaaic  coef- 
ficents  are  shown  in  Pl^.  3*19  to  3*22  for  an  aspect 
ratio  six  foil  and  in  Fig.  3*23  to  3.26  for  an  aspect 
rauio  one  foil. 

The  general  behavior  of  che  coefficents  with 
reduced  frequency  is  the  sane  for  the  two -dimensional 
foil  and  the  foils  of  finite  aspect  ratio.  The  magni¬ 
tudes  are  of  course  less  for  the  finite  aspect  ratio 
foils. 


REDUCED  FREQUENCY- k 


REDUCED  FREQUENCY- k 


FIG,  3.19  LIFT  COEFFICENT  DUE  TO  HEAVE  FOR  AN  ASPECT 
RATIO  6  SUPECAVITATING  FOIL  COMPARISON  WITH 
TWO-DIMENSIONAL  THEORY 
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CENTER  SECTION  Cm 
NUMERICAL  RESULTS 
TOTAL  Cm 

NUMERICAL  RESULTS 
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FIG.  3.20  MOMENT  COEFFICENT  DUE  TO  HEAVE  FOR  AN  ASPECT 
RATIO  6  SUPERCAViTATING  FOIL  COMPARISON  WITH 
TWO-DIMENSIONAL  THEORY 
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NUMERICAL  RESULTS 
TOTAL  Cl 
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REDUCED  FREQUENCY- k 


FIG  3,2!  LIFT  COEFFICENT  DUE  TO  PITCH  FOR  AN  ASPECT 
RATIO  6  SUPERCAVITATING  FOIL  COMPARISON 
WITH  TWO-DIMENSIONAL  THEORY 


REDUCED  FREOUENCY-k 


FIG.  3.22  MOMENT  COEFFICENT  DUE  TO  PITCH  FOR  AN 

ASPECT  RATIO  6  SUPERC.AVITATiNG  FOIL  COMPARISON 
WITH  TWO-DIMENSIONAL  THEORY 
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3,23  LIFT  COEFFICENT  DUE  TO  HEAVE  FOR  AN 

ASPECT  RATIO  ONE  SUPERCAV ITATI NG  FOIL 
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1.0  2.( 

REDUCED  FREQUENCY- k 


FIG,  3.24  MOMENT  COEFFICENT  DUE  TO  HEAVE  FOR  A 
ASPECT  RATIO  ONE  SUPERCAVlTATI NG  FOIL 
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REDUCED  FREQUENCY-  k 

FIG.  3.25  LIFT  COEFFICENT  DUE  TO  PITCH  FOR  AN 
ASPECT  RATIO  ONE  SUPERCAVITATIN6  FOIL 

1 19 


REDUCED  FREQUENCY- k 

FIG.3.26  MOMENT  COEFFICENT  DUE  TO  PITCH  FOR  AN 
ASPECT  RATIO  ONE  SUPERCAVITATING  FOIL 
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CHAPTER  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


Numerical  solution  of  the  linearized  Integral 
equation  for  a  sheet  of  pressure  slnguleucltles  has  been  shown 
to  constitute  a  powerful  tool  In  dealing  with  steady  and 
unsteady  problems  of  thin  lifting  surfaces  moving  through  a 
liquid.  The  general  technique  of  deriving  the  kernel  function 
analytically  and  solving  the  resulting  Integral  equations 
numerically  by  assumed  modes  results  In  an  efficient  procedure. 
It  is  felt  that  these  techniques  and  further  extensions  and 
refinements  constitute  the  best  linearized  solutions  available 
for  a  wide  varlty  of  lifting- surface  problems. 

Certain  comparison  with  expeixments  made  In  water  fall 

because  of  effects  which  are  outside  the  scope  of  an 

Invlscld  linearized  theory.  The  apparent  Inadequacy  of  the 

Kutta  condition  at  high  reduced  frequencies,  noted  In  Chapter  I, 

Is  an  example.  The  disagreement  between  experiments  and 

numerical  calculations  for  steady  siypercavltatlng  foils  at 

finite  angles  of  attack  could  also  be  mentioned.  As  yet,  there 

are  no  loading  measurements  on  oscillating  hydrofoils  of  large 

aspect  ratio  running  parallel  to  the  free  surface  at  low 
2 

values  of  kP  .  Experiments  on  supercavltatlng  foils  of  finite 
span  In  unsteady  vibration  are  currently  being  preformed  by 
Acosta  at  the  California  Institute  of  Technology. 
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The  method  for  fully  wetted  nonplanar  surfaces  has  wide 
applicability  to  Interference  problems.  Experiments  on 
steady  and  unsteady  motion  of  nonplanar  surfaces  such  as  T  tails 
would  be  very  valuable  to  check  the  accuracy  of  the  load 
predictions.  After  investigation  of  the  kernel  functions  for 
helical  motion,  this  approach  could  also  be  used  for  the  more 
complicated  problems  of  thin  lifting  propellor  and  helecopter 
blades.  The  treatment  of  the  effect  of  gravity  waves  on  an 
oscillating  foil  could  be  extended  to  compute  "gust  loads" 
on  a  foil  running  beneath  waves.  The  three’-dlmenslonal  kernel 
function  can  be  evaluated  numerically  without  difficulty  if 
finite-span  effects  are  of  Interest.  The  Influence  of  cavity 
closure  conditions  on  the  predicted  steady  and  unsteady  loads 
on  a  supercavltatlng  foil  should  be  investigated  further. 

Also  more  refined  models  for  supercavltatlng  foils  could  be 
developed  with  continuous  functions,  which  have  the  proper 
singularities,  assumed  for  the  pressure  source  modes.  The 
kernel  functions  for  a  supercavltatlng  foil  traveling  close 
to  a  free  surface  at  infinite  Proude  number  can  be  derived 
straightforwardly  by  means  of  images  and  the  numerical 
technique  extented  to  cover  this  case. 
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APPENDIX _ 1 


EVALUATION  OP  KERNEL  FUNCTK^S 
FOR  STEADY  AND  UNSTBAOT 
THREE  DIMENSIONAL  PLOV 


We  wish  to  evaluate  the  following  functionat 
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Evaluation  of  Ki 
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The  expression 
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is  a  generalized  function  defining  the  delta  function  except 
for  a  constant  of  proportionately » A  •  This  constant  is 
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We  define  two  functions  to  be  evaluated  separately 
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Following  the  method  of  ref* (6)  for  a  similar  kernel 
function  we  make  the  following  substitution  in 
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Preforming  the  indicated  differentiation  we  obtain 
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But  rroa  the  evaluation  of  we  have  found  that 
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The  evaluation  of  these  integrals  nay  be  found  is  Ref.  (13) 
(Bateman) . 
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We  make  the  following  definitions: 
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Evaluation  of  H. 


means  of  the  substitution  ^  =  y^ySthho  in  the  Integral 
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definition  for  H2  r  obtain  the  following  result: 
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We  recognize  as  a  generalized  function  defining  the 
delta  function  within  a  constant  of  proportionality.  To 
solve  for  this  constant, we  require  that 
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We  evaluate  the  integral 
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Appendix  11  The  Kernel  Function  for  Nonplanar  Surfacee 


The  kernel  function  for  a  nonplanar  surface  In  an 
Infinite  fluid  la 

I  S.  ^  ^  ^ 

The  correction  to  the  kernel  fimctlon  for  a  nonplanar 
surface  near  a  free  surface  at  Infinite  Proude  nunber  la 
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If  we  make 

the 

following  definitions 
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Then  the  kernel  function  In  equation  (ll.J^  becomes 


And  the  correction  terra  In  equation  (ll. 2)  becomes 
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<Fjj)  ■+ 

r/  <5/-  4’-  ?/'•-  (11.1?) 

The  functions  In  equations  (11.12)  to  (11.14)  can  be 
evaluated  using  ref.  (12)  and  ref.  (7),  They  are: 
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The  coablnatlon  of  these  functions  as  In  equation (ii.io) —  (ll.l4) 
will  give  the  very  complicated  form  of  the  kernel  f motion. 


The  behATior  of  the  kernel  function  near  the  line 
z  >  y  may  be  found  by  expanding  the  expression 


(11.19) 


for  saall  T,  r-  n  Is  the  normal  to  the  surface, 

t  Is  the  tangent  to  the  surface. 

For  small  r 


(11.25) 
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APPENDIX  111  Transformation  of  the  Fourier  Inversion  Integral 

We  consider  the  transformation  s«s  ,  ^  Ts  applied  to 
the  Fourier  Inversion  Integral 
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■nils  Integral  Is  also  written 
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Applying  the  transformation  s-s,  p-  Ts, 
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Let  s*-s  In  the  second  Integral  In  equation  (111.3) 
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¥e  then  let  'Y-  -T'  In  the  second  integral  of  (111,4) 


^"'1'’ iirf  (  sh  tlT- 


■i-r/ 

-  »  -  C» 


cSKri  ^ 


(111.5) 


Written  In  more  compact  form,  this  becomes 
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Iherefore  If, 
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then 
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